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A phenomenological prediction for radiative neutron capture is presented and compared to re¬ 
cent compilations of Maxwellian average cross sections and average radiative widths. The basic 
parameters for it, photon strength functions and nuclear level densities near the neutron separation 
energy are extracted from data without an ad-hoc assumption about axial symmetry - at variance to 
common usage. A satisfactory description is reached with an astonishingly small number of global 
fit parameters when theoretical predictions on triaxiality are inserted into conventional calculations 
of radiative neutron capture. These predictions (constrained HFB calculations with the Gogny D1S 
interaction) were tabulated recently for a large number of nuclei. For the photon strength a param¬ 
eterization of GDR shapes by the sum of three Lorentzians (TLO) is extrapolated to low energies. 
Level densities are influenced strongly by the significant collective enhancement based on the break¬ 
ing of shape symmetry. In the predictions for both the replacement of axial symmetry by the less 
stringent requirement of invariance against rotation by tt(R w ) leads to a global set of parameters, 
which allows to cover the range in nuclear mass number A from 50 to 250. The impact of non-GDR. 
modes adding to the low energy slope of photon strength is also discussed including recent data on 
photon scattering and other radiative processes. These are shown to be of minor importance for a 
comparison with experimental cross sections for neutron capture by even target nuclei. Here the 
triple Lorentzian (TLO) fit method for a parameterization of giant dipole resonances is normalized 
in accordance to the dipole sum rule and the droplet model with surface dissipation accounts well 
for positions and widths without additional parameters. Thus a reliable prediction for compound 
nuclear reactions also outside the valley of stability is expected. 


I. INTRODUCTION 

The ongoing discussion about triaxial shapes in 
heavy nuclei, recently often studied far off stability, may 
provoke the question, how well the widely used assump¬ 
tion about axial symmetry of most less exotic nuclei 
is founded on sufficiently sensitive experimental data. 
In the sense of formal logics conclusions from masses 
and level energies seem to be an unsufficient condition. 
The splitting of the Isovector Giant Dipole Resonance 
(IVGDR) is used as an indicator of axial deformation; 
a respective proof needs an adjustment of the apparent 
width locally for the five isotopes regarded 0], and a more 
rigorous approach appears helpful. In the present work 
this feature will be reviewed for 23 nuclei (as examples) 
in a wide range of mass number A under the assumption 
of broken axial symmetry. It enables a restriction to only 
one global fit parameter, which is valid for all nuclei re¬ 
garded, to parametrize the IVGDR-width. As the energy 
integrated absorption cross section is fixed by a sum rule, 
no fit of it is needed. A similar reduction of the number of 
free parameters is found for the other topic of this paper, 
the accordance of capture resonance spacings to a Fermi 
gas prediction. Again the assumption of broken axial 
shape symmetry strongly modifies the analysis of the ex¬ 
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perimental information as obtained from neutron capture 
by 140 spin-0 target nuclei. As an introduction to the 
question of nuclear shape symmetry, information from 
low energy nuclear structure studies in support of triaxi¬ 
ality will be presented first. Detailed Coulomb excitation 
studies have indicated the importance of the breaking of 
nuclear shape symmetries i.7] in accord to predictions of 
electromagnetic nuclear properties as clearly influenced 
by this breaking fg[. These studies go beyond the obser¬ 
vation of energy spectra, which have been interpreted in 
the past preferentially assuming at least axial symmetry 
of the nuclear shape 01M3. And they are superior 
to the electromagnetic excitation of only the lowest 2 + - 
level as a sufficient measure of axial deformation El- 
Interesting insight is gained by a comparison of various 
experimental data to recent theoretical work [l| allowing 
for the replacement of an 'adhoc' requirement of axial- 
ity by the less stringent assumption of invariance against 
rotation by an angle 7r, R k . 

The paper will first review the various connections be¬ 
tween nuclear shapes and electric quadrupole moments 
and transition rates; in the past these have significantly 
contributed our knowledge about fundamental nuclear 
properties. After a short discussion of electromagnetic 
sum rules for the nuclear dipole strength, a parameter¬ 
ization of the IVGDR as seen in photon reactions with 
heavy nuclei will be presented. In contrast to previous 
work the electric dipole strength is derived from IVGDR 
data without an assumption about spherical or axial sym¬ 
metry as usually made. The results will be regarded with 
respect to their impact on the radiative capture process, 
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for which the low energy behavior of this strength is of 
major importance. As was noted recently for nuclei with 
mass number A > 70 0 , the low energy tail of the elec¬ 
tromagnetic dipole strength function depends sensitively 
on the proper description of the shape of the IVGDR, in¬ 
cluding not only its position, but in particular its width 
and eventual split. On the basis of calculated [lj deforma¬ 
tion and triaxiality very few global parameters suffice for 
a parameterization of the dipole strength of heavy nuclei. 
The theoretical shape parameters will be used to evalu¬ 
ate the electromagnetic strength as derived from pho- 
tonuclear data in the IVGDR region. Photon emission 
and absorption data are important here fl4j - fl9j as they 
also deliver information on the possible influence of other 
multipolarities. The breaking of nuclear shape symmetry 
on the level density will also be regarded, as it strongly 
influences the phase space for radiative neutron capture 
and hence the size of the cross section for this impor¬ 
tant process. The radiative capture of fast neutrons by 
heavier nuclei plays an important role in considerations 
for advanced nuclear systems 011 and devices aiming 
for the transmutation of radioactive nuclear waste. This 
process is of interest also for the cosmic nucleosynthesis, 
especially for scenarios with such high fluxes of neutrons, 
that their capture becomes the main production process 
to reach heavier nuclides beyond Fe j23|, 0. As the ex¬ 
perimental studies which are the basis for respective pre¬ 
dictions can mainly be performed for nuclei in or close 
to the valley of beta-stability, a small number of global 
parameters increases the reliability of any extrapolation. 


II. NUCLEAR QUADRUPOLE MOMENTS AND 
DEFORMATION 


The electromagnetic response of nuclei has played an 
important role for the exploration of the size of nuclei 
and the departure of their shape from spherical symme¬ 
try was first indicated by a splitting of atomic transitions 
due to the nuclear electromagnetic field [25]. Much im¬ 
proved and accurate hyperfine structure measurements, 
partly using laser techniques, determined the ‘spectro¬ 
scopic’ electric quadrupole (A = 2) moment Q s of the 
ground state 0 in nearly 800 odd nuclei. In addi¬ 
tion, Q s -values from the reorientation effect in Coulomb 
excitation as well as from muonic X-ray data became 
available for many even and odd nuclei. In even nu¬ 
clei, Coulomb excitation-reorientation data can also yield 
‘spectroscopic’ quadrupole moments for excited 2 + -states 
Q s (2 + ), as compiled [26| for nearly 200 isotopes. In many 
cases the sizes and signs - positive and negative values 
are observed - of these quadrupole moments were not 
in agreement to predictions for single particle or hole 
configurations @ and a seemingly obvious picture [9j 
interprets this as the result of the rotation of a non- 
spherical body with an ‘intrinsic’ quadrupole moment 
Qo = —7/2Q s (2 + ). If the sign of Q s (2 + ) was determined, 
it mostly was negative indicating an oblate charge dis¬ 


tribution in the laboratory system and an intrinsic pro¬ 
late deformation. For most even nuclei away from closed 
shells the large Qo values observed were interpreted as a 
strong indication for a cigar like deformation, especially 
for lanthanide (A « 170) and actinide nuclei (A ss 240). 
Mainly for isotopes in the near magic Os-Pt-Hg region 
positive Q s (2 + ) were observed 0, indicating oblate in¬ 
trinsic shapes, i.e. Qo < 0. Assuming axial symmetry 
and a homogeneous distribution of the charge within the 
nuclear volume, a rotational model was formulated, in 
which the intrinsic electric quadrupole moment Qo of 
even nuclei is related to the radius difference A R be¬ 
tween the long and the two short axes of the shape by 

i: 


Qo = \j^-ZR 2 0{\ + b0)\ 
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The rotational model relation Q] between deformation 
f3 , Qo and A R is widely applied when electromagnetic 
data are related to calculated nuclear (mass) deforma¬ 
tions usually characterized by /?. For years b ss 0.16 
was used [3, 0, but reference is also often made to a 
compilation of electric quadrupole transition widths 0, 
which proposed b = 0 as approximation. In the same 
paper b = 0.36 was applied in a more detailed compar¬ 
ison to calculations, performed assuming zero hexade- 
capole moments. Besides this ambiguity in b several def¬ 
initions proposed as deformation parameters in the litera¬ 
ture 0,1 iMSI differ from the ‘standard’ definition as 
given here. Not all of them - even the ones in widespread 
use - conserve the volume when departing from spheri¬ 
cal symmetry, especially for large /3 or Qo- This may 
lead to unwanted and often neglected contributions from 
compressional energy in respective calculations. 

Intense experimental investigations .§, [H, 0 concern¬ 
ing the energies and angular momenta of the levels in 
heavy nuclei have been regarded in view of single parti¬ 
cle as well as collective excitations like rotations or vi¬ 
brations of the nuclear body as a whole. In the data 
analysis and in theoretical studies (6] spherical or axial 
symmetry was usually assumed, albeit a sensitive exper¬ 
imental verification was missing. Here observed level en¬ 
ergies and differences between them, eminent in gamma 
transition energies, have played an important role. But 
eventual strong energy shifts due to configuration mix¬ 
ing may confuse the arguments making the enhancement 
of electromagnetic strength a more reliable indicator of 
nuclear collective motion. Subsequent to a more formal 
discussion of the influence of nuclear shapes on photon 
interactions data will be presented indicating how break¬ 
ing axial symmetry influences electromagnetic strength 
and photonuclear reactions. Finally it will be shown how 
resulting new degrees of freedom enhance level densities 
and consequently radiative capture cross sections. 
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III. QUADRUPOLE TRANSITION 
PROBABILITIES AND TRIAXIALITY 

The enhancement seen in experimental data on electric 
quadrupole (E2) transitions from the ground state 0 
over predictions for a transition to a configuration formed 
by exciting only one or a few single particle orbitals is 
linked 0 to the breaking of spherical symmetry in heavy 
nuclei away from magic shells. The model of a rotating 
axially symmetric liquid drop with a quadrupole moment, 
representing an even nucleus, predicts one 2 + -state with 
a “collective” i.e. enhanced E2-transition width. In this 
model the intrinsic structure for the ground state 0 + and 
the lowest 2 + -state r are assumed to be the same. Then 
the reduced matrix element (in e fm 2 ) of the electric 
quadrupole transition E 1 (in MeV) of this ‘rotational’ 
state r to ground is related 0 : §] to Qo (in fm 2 ) by: 

|(r||E2||0)| 2 =B(E2,0^r) = 21_Q 2 (2) 

The E2 ground state decay width T r 0 (E 7 ) (in MeV) is 
obtained from the relation (valid in general): 

4tt a 

r rO (E 7 ;E2) = --^|(r||E2||0)| 2 ; (3) 

2 J r -\-1 
9 = 2 J 0 + 1 

where a ei h and c are the fine structure constant, Planck’s 
constant and the velocity of light; Jo and J r are the spins 
of the ground state and the excited level. The reduced 
transition probability B(E2) used in [2] is often expressed 
in e 2 fm 4 , but in this work we prefer to divide the E2 
operator directly by the elementary charge e, and do the 
equivalent for the B(EA). The Qo and B(E2) are observ¬ 
ables, and we will use them to express an experimentally 
determined deviation from sphericity, whenever possible; 
the deformation /3 is a model parameter, and in view of 
ambiguities in the definition of /? it will mainly be used in 
relation to theoretical calculations and as defined there. 
To describe the rotation related quadrupole transition 
between the ground and the lowest 2 + -state one often 
assumes ad hoc a uniform axial charge distribution; the 
same assumption has been made in the past 0 to relate 
axis ratios derived from IVGDR shapes to Qo or /?. The 
axis ratios used in this study are taken from the calcu¬ 
lations flj regarded as reliable source for the amount of 
triaxiality. 

One serious shortcoming of the axial rigid rotor model 
is the fact, that it only predicts one ‘collective’ 2 + - 
state. Experimentally at least two 2 + -levels with en¬ 
hanced transitions to the ground state are observed in 
nearly all even nuclei. This has led to the assumption 
@, Ha| that the coupling of the collective rotation to a 
collective quadrupolar vibration around a deformed ba¬ 
sis state has to be invoked; it can either be along the 
symmetry axis (/3 -vibration) or perpendicular to it (7 
- vibration) 0. If for a respective Hamiltonian a so¬ 


lution is found, which describes the level scheme simul¬ 
taneously to the B(E2)-values reasonably well, one may 
consider this a determination of the mean values of /3 
and 7 , but for their dynamic uncertainty additional con¬ 
siderations are needed. Reference is often made to more 
microscopic calculations in a self-consistent scheme and 
various approximations [2, |29M32| have been proposed re¬ 
cently. The first of these calculations is available for prac¬ 
tically all heavy even nuclei between the neutron and pro¬ 
ton drip lines. In Fig. |Tj the correlation between 7 and 
Qo is depicted for nuclei in the valley of stability as it 
results from this CHFB-calculation. Assuming only/Q- 
invariance they find non-zero triaxiality ( 7 ) 7 ^ 0 for many 
nuclei, and in some cases the predicted standard devia¬ 
tion does not include 7 = 0, i.e. cos( 37 ) = 1 . In Fig. 
[I] the density of symbols is of significance, as all the nu¬ 
clei from a small band near beta-stability are depicted. 
The clear clustering at Qo ~ 200fm 2 and cos( 37 ) ~ 0.2 
will play an important role for the discussion of IVGDR 
data in these numerous nuclei, previously often regarded 
transitional between axial and spherical in shape. 



FIG. 1. (Color online) Correlation between cos( 37 ) and Qo 
in « 170 even nuclei with 60 < A < 240; the respective data 
are taken from a CHFB-I-GCM calculation !l| for nuclei in 
the valley of stability and correspond to its minimum and 
to nuclei between two less and two more neutrons. The bar 
lengths represent the standard deviations in 7 as given by 
these calculations and tabulated as supplemental material; 
the dashed blue curve serves as eye guide to compare the 
trend of the calculations to experimental data as presented in 
Fig. El where it is depicted as well. 

The CHFB calculations used here “are free of param¬ 
eters beyond those contained in the Gogny D1S inter¬ 
action” (adjusted to the properties of nuclei with small 
Qo) and based on “a density-dependent HFB approxima¬ 
tion. They describe simultaneously the gross properties 
depending on the average field as well as the fine effects of 
the pairing correlations via the Bogolyubov field with the 
same force” [lj. The use of constrained wave functions 
allows to project on good angular momentum by the gen- 
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erator coordinate method, which is important at low an¬ 
gular momentum. Information on triaxiality obtained for 
the intrinsic system in Hartree-Fock-Bogolyubov (HFB) 
calculations is subject to significant change when the or¬ 
der of variation and projection on angular momentum 
in the observer’s frame are interchanged Q], fddl ]. Long 
ago it was pointed out [HJ, that the quantum mechani¬ 
cally proper variation after projection may shift the 7 - 
oscillation centered at axiality to ( 7 } 7 ^ 0. Another 
way to circumvent this problem is to parameterize the 
quadrupole degrees of freedom as seen by the observer 
in the laboratory, as is done implicitly by a interacting 
boson approximation (IBA) [35l - l37| . A rather convincing 
agreement to E2-transition data as well as level energies 
is achieved by this group theoretical ansatz when adjust¬ 
ing parameters for a given region of the nuclide chart and 
then comparing experiments to the predictions for other 
nuclei in that region. A global description of all heavy 
nuclei was not reached yet, but by a distinction between 
neutron and proton ‘Bosonic’ modes an account for tri¬ 
axiality was made (IBA-2) {38]; an inclusion of Giant 
Resonance modes would be of great interest. Another 
way to circumvent the problems related to the necessary 
projection into the observer’s frame is the construction 
of rotation invariants from experimentally observed tran¬ 
sition rates j§). Data from of heavy ion induced multiple 
Coulomb excitation of low lying levels in heavy nuclei 
were analyzed on the basis of such invariants [7j, [39( and 
for many heavy nuclei the breaking of axial symmetry 
was indicated by these experimental data. This analysis 
excels an older, comparatively simple, model of a rigid 
triaxial rotor [40| which directly delivers two collective 
2+-levels. Generalizing Eq. ©) Q 0 is replaced by Q,; de¬ 
fined by the sum of several (in the case of a rigid triaxial 
rotor 2) squared E2-matrix elements: 

£ I ( r ll E2 ll°> I 2 (4) 

r=l,s 

Experimental data |4lj| show that a limitation to one term 
in the sum (s = 1 ) leads to an error of less than 10 % in Qi 
and thus for the figures it suffices to use Eq. © instead 
of Eq. © to depict Q derived from B(E2)-values. In 
the rotation invariant ansatz Q the deviation from axial 
symmetry is described by the parameter cos( 37 ), which is 
in principle directly related to transition rates observable 
in multiple Coulomb excitation saiiffl. In quantal 
systems like nuclei only expectation values are accessible 
to measurements; Qo and cos( 37 ) in Figs.© [2]) are to be 
understood as symbols for the actually determined quan¬ 
tities; explicitely these are {Qi) and (Qf cos( 37 )) . In 
various nuclei the inclusion of triaxiality produced tran¬ 
sition rates close to data [~T~iT - TfGl | , especially when surface 
vibrations are included as a perturbation. But in most 
heavy even nuclei more than 2 collective 2 + -levels are 
observed and this scheme has to be generalized further. 
Employing tensor algebra and assuming reflection sym¬ 
metry as well as identical distributions of protons and 
neutrons Eq© was shown Bam to be valid also for 


heavy nuclei in general, as long as the sum includes all 
relevant 2 + -strength, especially all collective 2 + -levels. 
Then Eq©) with s > 2 represents an un-weighted sum 
rule which is rotation invariant and thus valid in the lab¬ 
oratory as well as in the intrinsic frame. 

In Fig. [ 2 ] the correlation of the parameters Qo and 
cos( 37 ) is displayed for more than 150 nuclei, for which 
relevant data have been published (1. B. liffiL l39l . l47l— 

I 49 I ]. The rather complex experimental investigations on 
triaxiality were performed for a limited number of nu¬ 
clei only, but a comparison to the prediction jl| shown in 
Fig. ©indicates a clear trend to triaxiality with decreas¬ 
ing Qo is obvious. Here (cos( 37 ) —> 0), whereas most 
well deformed nuclei and especially the actinides show 
smaller deviation from axiality. The small number of nu¬ 
clei which are oblate already at low E x (cos( 37 ) < 0) 
does not allow similar conclusions, and for very small Qo 
a definition of triaxiality is not appropriate. 



-4 -2 0 2 4 6 8 10 12 

Q 0 /100 (fm 2 ) 


FIG. 2. (Color online) Correlation between cos(37) and Qo 
in ~ 150 even nuclei with A > 60, for which respective ex¬ 
perimental data are available. The bars correspond to ex¬ 
perimental uncertainties and the dashed blue curve serves as 
eye guide to compare the trend of the data to the one of the 
calculations as presented in Fig. © where it is depicted as 
well. 

If not determined directly, the sign of the latter was 
assumed to coincide with the one of cos( 37 ); unambigu¬ 
ously opposite signs for these two quantities were not 
reported yet. Data for odd nuclei [50] are omitted in the 
Figure, we only state that they support the findings pre¬ 
sented here. Also not given in the Figure are variances 
that can be derived in principle from an extension of the 
above formulae [4, jj, Hi HU , but the experimental uncer¬ 
tainties of respective data do allow a rough extraction 
of such information only. The experimental data show a 
very similar trend of the axiality increasing with Qo as 
seen in the calculation. The trend as indicated as blue 
dashed curve in figs. [H and 0 suggests a speculation of nu¬ 
clear shapes represented by one parameter Qi only, with 
the axiality depending on it. 
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The split of the giant dipole resonances 0 is related 
to the axis lengths i?i, 2 , 3 - Assuming R n -invariance only, 
the combination of Eq. (13) to Eqs. (19-21) of ref. [ 8 } 
results in relations for the three axes Ri of the ellipsoid 
in the body fixed system as derived from its shape pa¬ 
rameters: 


5 Qi cos( 7 ) = Z(2Rl -R{- R%) 

5 Qi sin( 7 ) = V3Z{Rj - R%) (5) 

Rq = R1R2R3 

The last line uses the concept of the conservation 
of the nuclear density in an equivalent ellipsoid, which 
has the same volume V = 4/37t1?q, and the three har¬ 
monic oscillator constants w*, are inversely proportional 
to Ri. Assuming the same charge Ze for the equiva¬ 
lent ellipsoid as for the non-spherical nucleus, the Eqs. 
© can be related directly to the information about R 0 , 
the quadrupole moment Qi and the triaxiality 7 of that 
nucleus. The similarity between observations and the 
CHFB-calculations Jlj, as shown in Figs. |T]and[ 2 j sug¬ 
gests to use th ese as a reference. The equivalent sphere 
radius R 0 = -y/(5/3)(i? c ) with the charge radius R c , tab¬ 
ulated as well [l( , will also enter into the predictions dis¬ 
cussed below. In ref. [l] the tabulated /3 and 7 values 
are related to the oscillator parameter ratios Pchfb and 
Qchfb , the direct outcome of these calculations (with 
the latter not describing a quadrupole moment). These 
are used in the following Eq. © as axis ratios to obtain 
the radii Ri from Rq'. 


Ri = Rq(PchfbQchfb )^ 1 ^ 3 

R2 = PchfbRi ( 6 ) 

R3 = QchfbRi 


with 


x = 


2/3 + 1 

Pchfb = exp (-arv/ilsinq) 


Qchfb = exp x 


3 yfi . 

- cos 7 — — sin 7 


In ref. jlj] supplemental material is given for 1712 even- 
even nuclei including their deformation parameters, from 
which the corresponding ‘triaxial oscillator parameters’ 
can be derived, which are inversely proportional to the 
axis lengths. The extraction of these axis ratios differs 
formally from Eq. © as well as from the Hill-Wheeler 
formula used by us before [l3|. It can easily been shown 
numerically, that for the small Qi in the range of interest 
the differences are below 25% and hence not significant. 
For the calculation of fission barrier heights more exact 
prescriptions may be needed, but, as was pointed out 
already [5l], H2|, the consideration of triaxiality is even 
more important in that case. 


IV. PHOTON ABSORPTION SUM RULES AND 
GIANT RESONANCES 


The non-resonant interaction between photons and ob¬ 
jects of charge Ze and mass M is quantified by the Thom¬ 
son scattering cross section 8n (Z 2 a e hc) 2 /3(Mc 2 ) 2 with 
a e denoting the fine structure constant, c the velocity 
of light and h the Planck constant divided by 27r; for 
208 Pb it amounts to 0.02 fm 2 only. In addition to this di¬ 
rect process a photon of sufficiently high energy E y = E 
excites nuclei from the ground state resonantly; this is 
described by a Lorentzian centered at the resonance at 
E r with total width T r : 

~ T 2 ^T r 

r % (E r 2 -E 2 ) 2 + E 2 r 2 r [ ) 


The integral of the absorption cross section er 7 over the 
resonance, which has spin J r is denoted by I r o: 


Ao — 


a 1 {E)dE 


g(irhc) 2 T r0 
El ; 


2J r 1 
2 Jo + 1 


( 8 ) 


The width r r o is the partial width of the transition be¬ 
tween the resonant level ( E r , J r ) and the nuclear ground 
state (0, Jo). As described by Eq. © for multipole order 
A = 2, it is directly proportional to the electromagnetic 
transition matrix element; a respective relation exists for 

A = 1: 

r r 0 (Ey,El) = ^^^1 <r||El||0) | 2 (9) 


As derived from very general conditions like causality, 
analyticity and dispersion relations the interaction of 
short wavelength photons with nuclei of mass number 
A = Z + N can be ‘integrated up to the meson thresh¬ 
old’ analytically, leading to the energy-weighted sum rule 
of Gell-Mann, Goldberger and Thirring (GGT) (53j|: 


Iei = 


r 

Jo 


^ 27U 


= 5.97 


m„ 
ZN 

~A~ 


r (E)dE 

ZN 
~A + 
A 

+ 10 


A 
10 

MeVfm 2 


( 10 ) 


Here m n and m n stand for the mass of nucleon and pion, 
respectively. The first term in the sum is the “classical 
(TRK) sum rule” for electric dipole radiation [13] and 
the second “contains all of the mesonic effects” and is 
assumed [55| to be accurate within 30%. It was approxi¬ 
mated by assuming “that a photon of extremely large en¬ 
ergy interacts with the nucleus as a system of free nucle¬ 
ons” , and a correlation to hadronic shadowing was inves¬ 
tigated to be weak [55]] . Eq. (flOl) includes all multipole 
modes of photon absorption and no arguments fbil about 
the nuclear absorption of photons with energies above 
m^c 2 are needed. Absorption by the nucleons does not 
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contribute below £7, = m^c 2 , but nucleon pairs and espe¬ 
cially p-n-pairs are strongly dissociated by photons with 
20 < Ej < 200 MeV. The respective “quasi-deuteron ef¬ 
fect” has been derived from the expression valid for the 
free deuteron by correcting for Pauli blocking (57j . 

Photo-neutron data are available j58l| for 208 Pb up to 
energies above m,c , they are shown in Fig. [3] and com¬ 
pared on an absolute scale to respective theoretical ex¬ 
pressions: 

1. a Lorentzian like in Eq. 0 with parameters from 
Eq. (fT21) (width r r = 3.26 MeV and pole energy 
E r = Eq = 13.6 MeV). I r is normalized such that 
it’s integral agrees to the first term in Eq. m, 
and 

2. the expression for the absorption corresponding to 
the quasi-deuteron mode [571 ]. 

Apparently the cross section above 40 MeV is well re¬ 
produced on absolute scale and it arises that its integral 
is close to the second term in Eq. m- It is worth 
mentioning that a similarly situation was published [57:] 
for other nuclei. In contrast to elastic electron scattering 



FIG. 3. (Color online)Cross section of photo-neutron pro¬ 
duction data [58] on 208 Pb in comparison to a Lorentzian 
for the isovector IVGDR (full line, magenta) and the quasi- 
deuteron effect (blue dotted line). In 208 Pb a deformation 
induced widening can be neglected as will become obvious in 
Fig.m 

mm the cross section shown in Fig. [3] shows only weak 
signs of giant resonances of multipolarity other than El. 
The enhancement between 20 and 30 MeV observed to 
exceed the sum of the two curves may be identified with 
the contributions from T<-component of the IVGDR [dll ] 
and the IVGQR j5S, IfiOj- For the quadrupole giant reso¬ 
nances the energies are known in several nuclei from ex¬ 
periment with electron mm or hadron beams IgE [63| 
such that the respective sum rule G5 can be used to de¬ 
termine an upper limit in Eq. (11011 . In a recent paper on 
the Ml-strength in heavy nuclei measurements for many 


nuclei are reviewed 0 and from the quoted values for 
Ml-strength the contribution of it to Eq. ca will be 
discussed later. 

The photo-disintegration of nuclei is one of the nuclear 
reactions studied first ever; soon after it has been recog¬ 
nized as dominating photon absorption by an excitation 
of the IVGDR. This isovector nuclear electric dipole res¬ 
onance represents a strongly collective oscillation of neu¬ 
trons against protons. The good agreement in the region 
below 20 MeV seen in Fig. [3] indicates that the first term 
in Eq. m, a Lorentzian in accordance to the classical 
electric dipole sum rule (TRK, 0) is a good ansatz to 
describe the IVGDR. Although it is not a distinct quan¬ 
tum level, but a sum of densely packed levels resonantly 
enhanced near the ‘pole energy’ E r = Eq of Eq. 0 , it is 
described well. Derived from Eq. 0, but now applied to 
a wide giant ‘doorway’ forming an envelope over narrow 
electric dipole states excited by El radiation, Eq. m is 
used for the parameterization of the IVGDR. The pos¬ 
sibility of a Lorentzian was predicted theoretically 0] 
and tested numerically at hand of data for a few nuclei 
[65| . For the sum of k Lorentzians, with k characterizing 
an eventual split, the main term of Eq. m, the classi¬ 
cal sum rule, is divided equally into k fractions, assuring 
normalization: 


(Ufa 

dE 7 


(E 7 ) 


= G 


El,IV 
abs 


(Ey) 


= 5.97 


ZN 2 
A kir 


E 


(EI 



fm 2 


( 11 ) 


It will be shown in the next section, that the TRK sum 
rule is especially well fulfilled, when accounting for the 
breaking of axial symmetry. Especially the lanthanide 
and actinide nuclei appear to be quite close to axiality, 
as obvious in experimental data [6611671 ]. It will be shown 
in section^ to IVIII that the extraction of total dipole 
strength from data is not very sensitive details in the 
region of the maximum, as long as the sum rule is re¬ 
spected, and the same is true for the height of the low 
energy tail. In Fig. [3] here approximating 208 Pb as spher¬ 
ical even above 10 MeV,fc = 1 is used. For heavy nuclei 
only very few parameters for pole positions and resonance 
widths are required in a global fit: The two historic theo¬ 
retical treatments of the IVGDR predict its energy rather 
well for medium mass nuclei [68| , respectively for the very 
heavy ones By using concepts of the droplet model 
these two approaches were unified [70] and hence IVGDR 
centroid energies Eq{Z 1 A) are well predicted in the range 
50 < A < 254. The nuclear radius used here is derived 
from the mean charge radius ( R c ) taken from the CHFB 
calculations [I] as Rq = y/ (5/3) • ( R c ); symmetry energy 
J = 32.7 MeV and surface stiffness Q = 29.2 MeV are 
taken from the finite range droplet model 0- Only one 
additional parameter, an effective nucleon mass has to 
be adjusted [13| to m e ff = 800MeV/c 2 in an overall fit 
to the IVGDR data; due to the use of radii, deformation 
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and triaxiality from the CHFB calculations [lj it differs 
from our earlier work [13], where it was taken as 874 
MeV. By the use of Eq. © the energies of the three 
resonance poles are derived from the spherical centroid 
Eq and the well known proportionality between Ei and 
1/Ri. As previously [313] we thus use (with all energies 
expressed in MeV): 


he 

8 J 

A 2 

R 0 \ 

1 TO* 

4 NZ 


1 + u — e 


1 3 u 


e = 0.0768, m=(1 —£)-A" 1/3 
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With one parameter for the energies and one for the 
widths, both adjusted to be equal for all heavy nuclei 
with A > 50 [3 lil] one gets a good agreement to mea¬ 
sured resonance shapes, as will be shown in Figs. 6 to 25. 
In principle, the nature of the IVGDR does not directly 
allow for the determination of its width in analogy to 
Eq. ©, and the widening due to deformation has to be 
accounted for by using k > 1 in Eq. (ED- A strict distinc¬ 
tion between damping or spreading and the deformation 
induced splitting - to be discussed below - allows to ne¬ 
glect a photon energy dependence of the widths extracted 
from IVGDR data. Hydro-dynamical considerations [74j 
predict the dependence of the damping width T; of an 
IVGDR on its pole energy Ei to be proportional to El - 6 ; 
this exponent lies between theoretical values [75] for one- 
and two-body dissipation. Of course, the proportionality 
constant c w has an uncertainty related to the selection of 
nuclei, which are included in the fit. With the axis ratios 
available from the CHFB calculations an inclusion of all 
nuclides for which respective data have been tabulated 
[76} yields c w = 0.045(3). The width parameter is vary¬ 
ing with Ei and not with E 1 like in previous work jl3f . 
As the slope of a Lorentzian sufficiently far away from 
Eq is quasi proportional to Tj , its uncertainty enters in 
the radiative width nearly linearly. In our earlier work, 
a value of c w = 0.05 had been used [3, III; there single 
Lorentzians were adjusted to IVGDR data for 88 Sr and 
208 Pb, assumed to have one pole only. The new values 
for Eq and c w given here are based on the calculations 
P, 30f and also account for the influence of the GQR’s, 
to be regarded below. 

Individual local fits [3 [zl to the apparent widths 
seen in photo-neutron data (using k = 1 or 2 in Eq. (fill) , 
selected ad hoc to optimize the fit, result in a surpris¬ 
ingly large excess above the TRK sum rule. The result 
of such local fits as performed within the Reference Input 
Parameter Library (RIPL) project (3 [77} is depicted in 
the top part of Fig. 2] It differs somewhat from earlier 
data analysis @ which yielded an average of 105(8)% 
for 100 < A < 200 and a finding for six Sn-isotopes, 
where an overshoot of only 10(7)% was found |78j. The 
large unsystematic scatter reported from both local fits 


yield strong arguments against their use for nuclear as¬ 
trophysics, albeit proposed [79j in that field. 
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FIG. 4. (Color online) Panel (a) shows the energies (top) and 
widths (bottom, both vs. mass number) resulting from y 2 - 
fits to the IVGDR in heavy nuclei, as compiled recently [771 ]. 
The fits are based on one or two Lorentzians and two points 
per nucleus are shown, if a 2-pole fit led to a smaller y 2 . Our 
calculations with 3 poles as described in section IV and V 
(TLO) are depicted as drawn curves (in blue and magenta). 
In panel (b) the resulting GDR-integrals as obtained by the 
Lorentzian fits from ref. ]3l[ are depicted in comparison to the 
TRK sum rule (dotted line), which is surpassed considerably 
in most cases, whereas TLO obeys it by definition. 


The lower panel of Fig. [I] shows how the pole ener¬ 
gies and Lorentz widths from the local fits as published 
within RIPL-3 [3l[ (77) scatter as compared to the de¬ 
scription of the IVGDR shapes, allowing three poles and 
thus a ‘triple’ Lorentzian (TLO) parameterization. This 
is a preview to the predictions derived from the global 
fit procedure detailed in the next section. As seen from 
Fig. [4] the TLO-metlrod results in a smooth dependence 
on A which is modulated only due to variations in shape. 
And the comparison to various photon data in IV GDR 
peaks and their tails as presented in sections 0 toM 
shows that the new ansatz clearly agrees better with data 
than shown by previous works 66, [fi2l[ZZ|- But care is 
indicated with respect to a detailed comparison within 
narrow energy bins: For some nuclei near magic shells 



















the level-density is rather small up to the IVGDR peak 
(and far below a Fermi-gas prediction; cf. section ITXl) 
such that fluctuations are observed as long as the experi¬ 
mental averaging is comparatively small. Obviously, sev¬ 
eral components are responsible for the apparent width 
of the IVGDR in heavy nuclei: 

(a) Spreading into underlying complex configurations, 

(b) Nuclear shape induced splitting, 

(c) Fragmentation and 

(d) Particle escape. 

From calculations for heavy nuclei using the Rossendorf 
continuum shell model [80, |8J| the escape width (d) in 
the IVGDR region was shown to be clearly smaller than 
the spreading width (a) derived by the global TLO-fit; a 
good agreement to data shows that T in Eq. m depends 
on the pole energies only and one global fit parameter. 
For the concept of fragmentation (c) of the configura¬ 
tions belonging to e.g. the IVGDR a quantification of 
these configurations is needed. A detailed shell model 
calculation (H, @ f° r the nucleus 208 Pb is depicted in 
Fig. 0 which is based on a large number of configurations. 
From it one may conclude, that the parameterization of 
this IVGDR by a Lorentzian is justified as well as the 
width predicted as given in Eq. (fT2l) . It should be noted 
here, that this prediction, like TLO, is normalized to the 
’’classical” sum rule, the main component in Eq. m- 
The shape induced splitting (b) willl be regarded in the 
next section. 



E y (MeV) 

FIG. 5. (Color online) Dipole absorption from a shell model 
calculation for 208 Pb with configuration mixing j83| (black his¬ 
togram for El). For comparison the Lorentzian from Fig. [3] 
is superimposed as dotted curve (magenta); it corresponds to 
Eq. m reduced to one pole 

The agreement especially in the high energy slope as 
seen in Fig. [3] and Fig. [5] supports our concept of admit¬ 
ting no photon-energy-dependence of the width. Such a 


decrease was postulated by the KMF-model prediction 
for /ei allegedly derived from the theory for Fermi liq¬ 
uids [84j . A detailed relation to the work of Landau or 
Migdal is not given in that work, which explicitly states 
below its Eq. (29), that a ’’direct comparison” between 
the ’’spreading width” in a Lorentzian for the GDR and 
the one of the theory of Fermi liquids is ’’difficult... and 
not clear” [s3| . and they further assume ad hoc without 
any additional arguments, that they coincide. In this 
context it should be mentioned, that fundamental theo¬ 
retical arguments have been used to show, that ’’Landau 
damping is not the appropriate process for descri bing the 
damping of the low-multipole giant resonances” j75[. 


V. ISOVECTOR GIANT DIPOLE RESONANCES 
IN HEAVY NUCLEI AND TRIAXIALITY 

The coupling of dipole and quadrupole degrees of 
freedom in heavy nuclei has been discussed thoroughly 
0 H3, HH-Hl} and calculations within the interacting Bo¬ 
son model (IBA) have obtained reasonable fits to exper¬ 
imental data for two chains of even isotopes [ 88 ] and 
some other nuclei f89[ dH) . In the dynamic collective 
model (DCM) [80 l9l| calculations on the basis of hydro- 
dynamical considerations were successfully compared to 
selected data [93, 93]. The importance of the breaking 
of axial symmetry for the GDR shapes was not discussed 
in the literature cited here and not in most of related 
work. In theoretical papers on a few nuclei triaxiality 
was mentioned with respect to the IVGDR @@ , but 
a full coverage of a wide range in the nuclide chart is 
missing. The parameterization presented earlier fl3j and 
specified in this section is much less ambiguous concern¬ 
ing the mode coupling, but it has one advantage as it in¬ 
corporates nuclear triaxiality explicitly by setting k = 3 
in Eq. m For the resulting ‘triple’ Lorentzian (TLO) 
description the resonance energy Eo is modulated by the 
ratios of the axis lengths Ri in the spirit of Eqs. (ED and 
0 (setting Ei/E 0 = R 0 /Ri). This direct incorporation of 
triaxiality makes TLO differ from previous attempts to 
obtain Lorentzian fits to photo-neutron data for a large 
number of heavy nuclei 0, [TH], HJ HI, H3, [z3, HU- In 
many nuclei, especially those of intermediate Q o, the 
local fits presented there may lead to a seemingly bet¬ 
ter agreement, but often they require quite unreasonably 
large values for the width of the IVGDR and for the inte¬ 
grated strength in comparison to sum rule predictions. If 
triaxiality is accounted for in addition to the quadrupole 
deformation a good description of IVGDR shapes can be 
obtained jl3| without treating the strength and width 
as free fit parameters. The deformation induced shift of 
the three axis lengths Ri versus the equivalent radius Ro 
is obtained from Eq. 0 , which is based on a CHFB 
calculation 1 , [30|. As outlined in section uni the last- 
mentioned work lists values for the quadrupole deforma¬ 
tion /3 as well as the corresponding 7 for a large number 
of even nuclei. Hence, a global prescription to predict 
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electric dipole strength can be derived from it, not lim¬ 
ited to using local information for single nuclei. 

When a parameterization of the electric dipole 
strength in nuclei with non-zero Qi is aimed for, the con¬ 
tribution (b) of the list in the previous section has to 
be treated sufficiently well. As shown in the top part of 
Fig. SI the TRK sum rule, Eq. (1101) . disagrees in many 
nuclei to the Lorentzian fits [63, lZZ| performed for the 
data of each nucleus independently without account for 
the possibility of broken axial symmetry. In these fits 
the width parameter was adjusted for each isotope sep¬ 
arately to fit the peak region and for A between 90 and 
150 an especially large discrepancy is observed as well as 
wide fluctuation with Z and A of this apparent breadth 
indicating a non-systematic variation which is difficult 
to conceive within the spreading concept. A similarly er¬ 
ratic dependence of the integrated IVGDR strength on Z 
and A was reported to result from this approach of 
fitting the photo-absorption data locally. In some cases 
the integrated cross section overshoots the classical sum 
rule given by Eq. (flT)l first term) by up to 100%. Ap¬ 
parently the two problems named are closely related, as 
the resonance integral is proportional to the product of 
height and width. As proposed previously [3 Isl, H7] , a 
solution for this problem is found by allowing axial sym¬ 
metry to be broken; this point will now be examined in 
further detail. As shown in section m accurate nuclear 
spectroscopic data suited to determine both deformation 
parameters are available only for a limited number of 
nuclei. The recent CHFB calculation |l| delivers prior 
information for Eq. (0 above, inserted to obtain the 
resonance energies in the sum of Lorentzian functions in 
Eq. (HU). This procedure leads to a significant split¬ 
ting into three equally strong IVGDR components which 
increases with deformation. As seen in Fig. [2 for many 
nuclei the splitting between the three components is com¬ 
parable in energy to their widths and thus not directly 
obvious from the data alone, especially in nuclides with 
Qo ~ 200 — 300fm 2 . These are not rare as the clustering 
depicted in Fig. |U shows, but the calculated triaxiality 
demands that all three axes are accounted for explicitly. 
This quite simple consideration explains the significant 
rise in apparent breadth as seen in Fig. [6] for the even 
Sm-isotopes with N = 86 to N = 92; 144 Sm was not 
included because of the uncertain cross section for the 
( 7 ,p)-reaction, and in 148 Sm, previously often regarded 
as spherical, the deformation induced split leads to a 
low tail, of importance for the envisaged extrapolation 
to even smaller energies. Special care is needed for nuclei 
near closed shells: The CHFB calculations do not fully 
account for the very deep mean field potential in such nu¬ 
clei and thus they produce too much of collectivity [3C|. 
Following what is said there, a reduction for nuclei only S 
nucleons away from a shell a factor for the /3-deformation 
[lj of 0.4 + 6 /20 is applied for S < 10. This expression 
is used for protons as well as neutrons and the larger of 
the two correction factors is taken; it results in a reduc¬ 
tion of the predicted [l[ /3-values by 40, 30, 20 and 10% 
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FIG. 6. (Color online) Photo-neutron production cross sec¬ 
tion 148 Sm [ 9 U (a) to 154 Sm (d) in comparison to the TLO 
sum of three Lorentzians (drawn curve) with Ei indicated as 
black bars. The dashed (purple) curves indicate the effect of 
shape sampling [§3, Is3, [99] • As outlined in the next section 
the data were renormalized by a factor 0.9 and the calcula¬ 
tions were folded with a Gauss ian t o simulate an experimental 
beam spread of a = 0.3 MeV [lOCj ]. Panels (b) and (c) depict 


the situation for 150 Sm and for 152 Sm, respectively. 


for the isotopes 148 Sm, 150 Sm, 152 Sm and 154 Sm and the 
corresponding agreement to the IVGDR data is shown in 
Fig. [G] It is worth mentioning here, that the increase of 
width Tj with Ei as predicted from Eq. m reduces the 
maximum in the upper of the two peaks, such that it is 
not two times higher than the lower, as expected from 
the near axial prolateness of 154 Sm. 

The supplemental material based on the CHFB- 
calculations [TJ does not only list the mean values for 
the deformation parameters as obtained via the GCM, 
it also gives their variances as resulting from quantum 
mechanical zero point oscillation. The respective Gaus¬ 
sian distributions obtained thus allow an instantaneous 
shape sampling (ISS) as shown earlier Q for isotope 
chains Mo [9^] and Nd [99j|. There the impact of ISS on 
the height of the low energy tail and thus on radiative 
capture was demonstrated to be negligible. Results for 
the Sm chain as shown in Fig. [6] indicate a minor influ¬ 
ence as well: the drawn black curves correspond to the 
TLO-prediction and the dashed purple curves stem from 
the analysis with TLO + ISS. The impact of this effect 
for other nuclei needs further study; we expect it to be 
small. The self-consistent extension to non-spherical nu¬ 
clei introduces no extra free parameters in addition to the 
global ones of the Gogny-force [H and an attempt of a 
fit to the IVGDR energies and widths succeeds with only 
the four parameters introduced in connection to Eqs. m 
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and m, two of which are known from LDM mass fits. 
In the evaluations presented in the next section the only 
local parameters for individual nuclei are the axis ratios 
calculated by CHFB [lj, and no dependence of the widths 
Tj (c/. Eq. (fl2l) on the photon energy E 7 was allowed. 
Invoking widths varying only with the pole energies Ei 
(and not with A and Z) opens the possibility for a global 
prediction of nuclear photon strength also for heavy ex¬ 
otic nuclei and an extension to energies below the neutron 
emission threshold S n . This has the potential of consis¬ 
tent predictions for radiative capture processes, where 
full satisf actio n is not reached with presently available 
methods [l0l| . But, as will also discussed below, other 
modes than the IVGDR contribute to photon absorp¬ 
tion and their influence will be quantified in sections I VIIll 
and [x|. 


VI. ELECTRIC DIPOLE STRENGTH IN THE 
IVGDR AND BELOW 


The interaction of photons with heavy nuclei at en¬ 
ergies above the neutron separation energy Sn is mostly 
resulting in reactions of type ( 7 ,xn), but for a determina¬ 
tion of the complete absorption cross section the emission 
of p, a’s or fission might also be important. As will be¬ 
come obvious the photon absorption in nuclei with low 
level density below the IVGDR pole may be fragmented 
into many small peaks, overlapping at higher excitation 
energy only. The radiative neutron capture discussed 
later may proceed via this domain and it is thus indicated 
to use average photon strengths f\ for the description of 
electromagnetic decay by multipolarity A. This average 
quantity, the strength function f\{Er 1 ) 1 was introduced 
fl5j ] for electromagnetic processes. The relation to the 
average photon absorption cross section and the integral 
I x bs over a gi yen energy interval around E 1 is: 






( 7 Thc) 2 g e ffEi; 


T x 

abs 


(■jrhc) 2 g e ffE% x 1 Ag 


(13) 

To use the strength functions f\{E 1 ) for excitation as 
well as decay processes and thus connect photon scatter¬ 
ing to radiative capture and photonuclear processes one 
has to suppose them to be direction independent (Axel- 
Brink hypothesis) [l5j] . Eqs. [H [9] and m may then be 
used to directly relate f\ to the electromagnetic decay 
widths of the resonant levels r in the integration interval 
A e '■ 


sity p r = 1/D r . Increasing (r r7 (.E 7 )) and the reduced 
matrix elements for the transition to the ground state 
results in an increase of fi\)(E But decays not to 
ground have to be corrected for [14[, as they reduce the 
observed gamma-intensities and for that Eq. m may 
be applied to obtain average branching ratios, which to¬ 
gether with pf play a role in the comparison to data. 
Such decays take place between levels which are both 
excited and there is no direct way to study them start¬ 
ing from target ground states. But the average quantity 
fx is usually [15| assumed to be insensitive to details of 
the nuclear excitations and hence one approximates col¬ 
lective electromagnetic transition strengths of energy E 7 
by /a(E 7 ) to be independent of the energies E r and Ef ; 
together with the notion of f\ being valid for excitation 
and decay, this assum ption is also part of the Axel-Brink 
hypothesis floi floi. In the present ansatz an eventual 
dependence of decays into a region of elevated nuclear 
temperature is not accounted for by a modified strength 
function, but rather by averaging over that region using 
the relevant level density (which may be characterized by 
an apparent temperature (cf. section HXl) . 

Several facts have to be regarded when using experi¬ 
mental photo-neutron or photon absorption data to test 
the TLO parameterization as strength-function predic¬ 
tion: 

1. Considerable discrepancies were reported for ex¬ 
periments performed e.g. at different laboratories 
pi - dOfij : in some cases energy calibrations may 
differ somewhat. 


2. Photo-neutron data were often obtained by using 
quasi-monochromatic photon beams with a rather 
wide energy distribution, which is incorporated by 
folding the calculations with a Gaussian of width 
< 7 = 0.3 MeV, a valu e not as large as some recently 
made guesses [M EoH. An apparent widening 
near the IVGDR peak is partly due to ISS as dis¬ 
cussed above. Also in the case of a bremsstrahlung 
distribution used as photon source quite some un¬ 
certainty may arise [1091 ]. 


3. In a number of cases the ( 7 ,p)-channel exhausts 
a significant portion of the photo-absorption cross 
section [63, |73, |9Z| • ^ was shown in an earlier paper 


97[, how the eventual influence of all open channels 


on the extraction of the absorption cross section 
from the existing data is tested by Hauser-Feshbach 
calculations. 


h (£7) = "a 


■ r 7 


A E ^ £^ a+1 


(r r7 (E 7 )) 

D r E% x+1 


(14) 


The integration in Eq. m runs from E 7 — AE/2 to E 1 + 
AE/2 and for the sum in Eq. (fill) all levels within this 
interval are included; the quantum-mechanical weight 
factor g e f / will be discussed in section IVIII The aver¬ 
age level spacing D r at the upper of the two levels con¬ 
nected by E 1 = E r — Ef is the inverse of the level den- 


4. At higher energy the competition by the ( 7 ,2n) 
channel becomes important and that requires in¬ 
volved subt ract ion procedures [93|, which may be 
questioned |l07l . llOq . 

5. Most of the targets used contain unwanted isotopes, 
and when some of them have a lower S n (like many 
odd isotones) the low energy yield has to be cor¬ 
rected j97j . 
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6. Below and above the pole of the IVGDR contri¬ 
butions from the giant quadrupole (GQR) modes 
cause effects of some significance. 

7. In heavy nuclei various excited states may ex¬ 
ist aside from giant resonances, which cause ob¬ 
servable photon strength also below the neutron 
threshold. It adds to the IVGDR tail not_ always 
smoothly, as Porter-Thomas fluctuations [1031 ] may 
randomly create strong peaks in spectra from a 
quasi-continuum of weakly populated levels. 

In the IVGDR region averaged experimental photo¬ 
dissociation data as obtained with quasi-monochromatic 
photons from positron annihilation have been compiled 
and are available [6(| H3, |76]. Such data do not exist for 
all stable isotopes and for some nuclei photo-absorption 
has been studied by an absorption technique, which may 
serve for a consistency check in spite of its systematic 
error due to the need to subtract t he strong atomic ab¬ 
sorption. As indicated repeatedly |lo4 Il08l Ill0l4ll3( | . 
cross sections obtained at Saclay should be reduced. The 
necessary reduction is probably related to difficulties in 
the a nalysis o f multi-hit events in the neutron det ecto r 
array |l08l . Ill2( | . In accor danc e to a precision study OH, 
confirmed by results [U Ill4l | from the radiation source 
ELBE, the photo-neutron data of that origin are hence 
multiplied in this work by 0.9, considered as suitable for 
various A. 

Items 2 to 4 influence the representation of the IVGDR 
peak region, but their effect on the tail a few widths T 
below the pole is insignificant - in contrast to items 5 
and 7. Photon strength data for energies below S n are 
even rarer fl3| . Similar as in many photo-neutron studies 
experiments with quasi-monochromatic beams are per¬ 
formed to gain information on absorption from photon 
scattering: In early e xper iments the tagging system at 
Urbana was used [14, 115] and an increasing number of 
measurements are done at the laser-photon backscatter- 
ing facilities HlyS (High Intensity y-Ray Source) at Duke 
University |ll6l Ill7| and AIST (National Institute of 
Advanced Industrial Science and Technology) in Japan 
jl06| . The photon intensities from these are not well 
known and bremsstrahlung data eventually serve for nor¬ 
malization. Results from scattering have to be corrected 
for branching to other than the ground state and in var¬ 
ious cases the bremsstrahlung continua also feed higher 
excited levels and their decay yield has to be subtracted. 
In the case of less complex spectra use can be made of 
data from the decay to well isolated levels EH, EH- 
More sophisticated schemes as developed e.g .at ELBE 
fsil 1971 Il20l Il2l| are based on statistical considerations 
already formulated some time ago El, El- Due to the 
fact, that the electromagnetic strength is responsible for 
the absorption as well as the emission of photons, an it¬ 
erat ive proced ure can lead to a self-consistent solution 
fstii Il22l . EH- Here, level spacings enter which often 
were taken from predictions and extrapolations eventu¬ 
ally questioned; in view of new results 124J on level den¬ 


sities as outlined in section IX]a reduction of 30% in the 
fx as compared to published values was applied in the 
case of ELBE-data and mentioned in the respective fig¬ 
ure captions. In principle, also photon yields observed 
after nuclear collisions can deliver strength-information, 
when it can be normalized via an ”external” fix-point. 
In the case of resonant neutron capture this is realized 
by ” using the absolute gamma-ray intens ities due to cap¬ 
tured thermal and resonance neutrons” }l25| . A similar 
normalization via neutron capture, now with respect to 
level densi ties, was attempted for data taken at the Oslo 
cyclotron jl26l - ll28| with 3 He and 2 H projectiles. The 
method of relying on neutron capture resonance spac¬ 
ings requires the knowledge of the energy and spin de¬ 
pendence of level distances; here theoretical assumptions 
have to be used. 

The subsequent compilation presents in figs. @ to [H 
photon strength function data and a comparison to the 
TLO parameterization with three IVGDR pole energies 
induced by the deformation. Here parameters are taken 
from the available CHFB-calculations El and globally 
determined fit values for c w in Eq. m and an effec¬ 
tive mass quantifying EH the centroid Eq. The (black) 
dashed lines represent the prediction of /ei thus derived 
with the resonance integral from Eq. HID ( in accordance 
to the TRK sum rule) equally divided among the three 
poles of TLO. In all these figures the three are indicated 
as black bars at the energy axis. It is worth mention¬ 
ing that the absolute heights in the low energy slope are 
nearly unchanged by the splits, albeit the apparent peak- 
height depends on it. As shown previously [831 ] and as will 
be further detailed in section la the effect of the photon 
strength on radiative neutron capture is strongest in the 
tail region below the neutron binding energy S n . Here up 
to eight additional ” minor” strength components may be 
of importance; they will be detailed in section IVIIII For 
a specification of their energy, strength and breadth we 
have used published data as listed with the figures. But 
the available information is by far not detailed enough 
to derive a systematic parameterization to evaluate their 
eventual influence on neutron capture with high preci¬ 
sion. This is demonstrated by the experimental data 
inserted in the figures; the result of our adjustment of 
the relevant parameters to derive an approximate agree¬ 
ment is depicted by a full (blue) line. Like for TLO only 
globally fitted quantities enter the calculations for these 
plots; they are compiled in Table E The plots start at 3 
MeV as below Thomson scattering by the nuclear charge 
surmounts the IVGDR tail and the zero at E 1 = 0 in Eq. 
m- It should be stressed here again, that the presented 
ansatz does not aim for a full theoretical understand¬ 
ing of the coupling between the IVGDR to quadrupole 
modes or other excitations, but only for phenomenolog¬ 
ical prediction of photon strengths, which is global and 
hence extendable to many nuclides. In view of the al¬ 
ready previously observed [15] ’’difficulty of accounting 
for the bump with the aid of a smoothly varying strength 
function”, we clearly distinguish between the Lorentzian 
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tail and ’’minor” strength. From the figures this appears 
to describe existing data at least as well as previous at¬ 
tempts like the ones based on the KMF model [13, [§5j . 
But it implies that for certain energy regions the strength 
published was observed only because it was significantly 
stronger than the one outside of them, which we assume 
to have been buried in the experimental background ( e.g. 
due to Compton scattering in the photon detector). 

The nuclei discussed in the following compilation were 
selected such that photo-neutron data and results for en¬ 
ergies below S n are available; still different A, Z and Qo 
are well represented. We deliver here a comprehensive 
collection of the data of interest, most of them were taken 
from an electronically accessible data base [76j, for which 
they had been extracted from original work. To improve 
the visibility of the data points, some of the excitation 
functions have been re-binned to around 0.6 MeV/bin, 
and points of no significance were suppressed, e.g. when 
their uncertainty is comparable to the value. The sub¬ 
sequent figures will demonstrate what features of photon 
strength information can be derived experimentally, and 
how well systematic trends become visible. For nuclides 
as light as 54 Fe channels competing to ( 7 ,n) cause difficul¬ 
ties for the extraction of the absorption cross section and 
hence the photon strength. In view of an increase of this 
problem with falling A the reason for not extending the 
study to lower A becomes obvious. Photo-dissociation 



FIG. 7. (Color online) Photon strengt h for 54 Fe predicted 
by TLO (dashed line in black). Data Il29tl (crosses ( 7 ,p); 
and full circles ( 7 , 11 ) nearly coincide between 15 and 18 MeV, 
whereas at higher energy the ( 7 ,p) channel dominates and 
their sum is not far - on absolute scale - from the dashed 
line in black (TLO) and the full line in blue, depicting the 
influence of minor strength. Photon strength data at low E 1 
(x in magenta) are from 56 Fe as obtained from 7 -decay after 
He-induced processes [l3dl . and the dotted line in red results 
from a one pole fit to the photon absorption cross section 
derived from the data [77). 

data for 54 Fe [l29l [ are shown in Fig. [7] to be not far from 
the TLO prediction in the IVGDR, when the ( 7 ,p) and 


( 7 ,n) channels are added up. The f\ derived [l30| | at 
lower E-y from 7 -decay in 56 Fe seem to agree reasonably 
well with the sum of TLO and considerable extra ‘pigmy’ 
strength near 7-8 MeV (c/. Table U). But it should be 
mentioned, that between 6 jnd 10 MeV two published 
data sets lie clearly below nm and above [126 ], and this 
is a typical example for problems arising for the kind of 
survey study as presented here. The observed rise below 
4 MeV does not correspond to any prediction considered 
in the present study, but it may eventually be similar to 
Ml strength recently predicted for some heavier nuclei 
[l32j| . For 1113 MeV, where the spin-flip Ml strength 
is expected, no data are available. In 78 Se a significant 



FIG. 8 . (Color online) Photon strength function for 78 Se: 
the TLO-prediction is depicted as dashed black line and a 
full line in blue represents the approximate influence of minor 
strength. Photoneutron results (black dots Il33il l are shown in 
comparison as well as photon scattering data at lowi5 7 ; these 
(x-symbols in magen ta) were obtained with bremsstrahlung 
from ELBE [83, 1.2.'! and reduced by 30% in view of a ’’new” 
level density ansatz, cf. section ED. 

increase over the extrapolated tail is observed from pho¬ 
ton scattering investigated at the radiation source ELBE 
for 7 < Ej < 10 MeV, although these data [13 were 
reduced by a factor of 0.7. This is a correction for the 
disagreement of the level density ansatz used in the data 
analysis [§3} with the one published recently [l24| , which 
will also be discussed in section m The IVGDR data 
were deriv ed fr om photo-neutron experiments performed 
at Saclay [l33j . 

Photon strength func tions have been published previ¬ 
ously for 88 Sr [rJ il24l |. as well as for 92-100 Mo [97]; a 
width independent of E 1 was used, a good agreement to 
the TLO prediction is observed and the integrated cross 
section obeys the TRK sum rule. The similarity between 
the low energy slopes in the experimental data of all Mo 
isotopes led to the suggestion [lj, [93] of an extrapolation 
with energy independent width under the assumption of 
broken axial symmetry. Looking at 98 Mo in the tail below 
10 MeV shows that bremsstrahlung data corrected sta- 
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FIG. 9. (Color onlin e) Photon strength for 98 Mo from pho¬ 
toneutron data Il34l1 (black circles) in comparison to the 
sum of three Lorentzians (TLO) as described for Fig. [5] . 
The data below 9 MeV are from elastic photon scat tering 
98 Mo(7, 7) observed with monochromatic photons [13511 (blue 
x) or bremsstrahlung tl20j ] (magenta + symbols, modified by 
0.7). 


tistically for inelastic scattering [971, Il20j are very close 
to ne wer results obtained with quasi-monochromatic pho¬ 
tons [l35| . At variance to this agreement photon strength 
function d erive d from gamma decay after inelastic 3 He 
scattering |l27l | are nearly a factor of 3 smaller at 7 MeV 
as compared to the data as shown in Fig. [9] and this 
may indicate, that the so-called Oslo method may suffer 
from the uncertainty in level dens ity mentioned above. 
The recently revised data fl0fii . H36j agree around 7 MeV 
with the photon scattering results and thus also with the 
TLO prediction with minor strength added. These new 
data overshoot at higher E 7 and are lower by ss 40% 
near 3 MeV; similar to Fe an increase of strength was 
observed at low energy and shown to be in accord to a 
shell mo del c alculation for Ml-strength in the nuclides 
94 “ 96 Mo [l32j| . New photo-neutron data below 11 MeV 
[l06l ] with smaller error b ars a nd smooth dependence on 
Ej replace the older ones [mj. 

Quasi-elastic photon scattering from natural Sn has 
been studied long ago EH at the tagging set up in¬ 
stalled at Urbana and ’’intermediate structure” in ad¬ 
dition to the IVGDR tail has been identified. The ab¬ 
sorption cross sections were derived from scattering data 
[l4[; their branching correction by inserting constant av¬ 
erage resonance widths may overestimate a a bs by at most 
20% EH- To improve the overlap with the ( 7 ,n)-data a 
reduction by 0.8 was applied in Fig. [TO] and the result¬ 
ing values are shown together with the cross section for 
118 Sn ( 7 ,xn) [l37j . obtained with positron annihilation in 
flight; the surprisingly large strength near 10 MeV may 
be related to a target admixture of odd isotopes, simi¬ 
lar to what was worked out for Mo targets [97]. Recent 
experiments with laser backscattered photons 1061 ] sup¬ 


port such an assumption. From a high resolution photon 



FIG. 10. (Color online) Photon strength for 118 Sn from 
photon scattering fl4| (blue x symbols ) at low and photo¬ 
dissociation [1371 ] (black circles) at higher energy in the 
IVGDR range. The data are compared to TLO (see the cap¬ 
tion of Fig. [5]). Recently reanalyzed data [ 5 } from gamma 
decay after 3 He induced reactions are shown as + symbols in 
magenta with large uncertainty bars. 

scattering experiment EH with correction for branching 
losses a strength enhancement near 6.5 and 8 MeV was 
reported for 116 Sn and 124 Sn, similar to what is shown 
here for 118 Sn. A recent study of 112 Sn and 120 Sn |138| 
uses statistic al co rrections for inelastic scattering as pro¬ 
posed earlier [l2C)t | and again finds similar dipole strength 
as reported earlier. In addition the last-mentioned inves¬ 
tigation has detected by a fluctuation analysis the need 
to increase the final photon strength by nearly a factor 
of two with respect to the sum of peaks observed in the 
spectra and resolved with a: 2 keV (FWHM) energy res¬ 
olution. Whereas data on / 7 derived at Oslo for various 
Sn isotopes and other nuclides with 3 He-induced reac¬ 
tions show significant differences from the photon scat¬ 
tering data shown here, the recent reanalysis for 118 Sn 
[l39l ] does not - as depicted in Fig. [TU] The overshoot 
seen at the higher energies can be explained by the in¬ 
fluence of the giant quadrupole resonance and is of no 
importance for radiative neutron capture. 

In Fig. HU results from photo-neutron emission from 
" at Te, multiplied by 0.9 - as done in general for data 
from Saclay - agree well to TLO above the IVGDR peak. 
The various isotopes in the target may cause a widening 
of it and below the isoscalar component of the GQR is 
expected from the systematics for this mode; its influ¬ 
ence on the photoneutron cross section is not completely 
clear. The low energy data support the finding of extra 
strength near CI.42.Eo, also seen in Figs. [TU] and [12] and 
to some extent in the data below 10 MeV for all nuclei 
presented here. This ” inter med iate structure” strength 
was observed since long Ei EH by photon scattering in 
many heavy nuclei - and labeled pigmy strength. It may 
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partially overlap Ml-strength 53 seen in various nuclei. 
Such peaked extra photon contributes little to the dipole 
sum in Eq. m at variance to a decrease of Trv n dr wit h 
decreasing E 1 , as favoured for long time [73, S, Il25j | . 
Whereas this KMF-ansatz often disagrees to the dipole 
sum rule, the TLO-approach presented here always obeys 
it and still agrees to measurements. For 138 Ba scattering 



FIG. 11. (Color online) Photon strength for 130 Te calculated 
as sum of three IVGDR-Lorentzians (TLO)(see the caption of 
Fig. |8]) in com paris on to data from photoneutron production 
in natural Te Il37ll (black circles). The data below S n (x, 
magenta) are from a careful analysis ll40H of scattering data 
taken with quasi-monoenergetic photons, yielding lower and 
upper limits for the photo-absorption cross section. 

data have been taken with quasi mono-e nerge tic photons 
at the laser backscattering beam at HlyS [l!7j , whereas a 
bremss trah lung experiment at ELBE was performed with 
136 Ba 5H1 ■ A s one can assume that the two data sets 
result in very similar absorption cross sections, they are 
show n tog ether in Fig. [12] A possible increase below 3 
MeV 122] may be of similar origin as obs erva t ions in Fe 
and for N ss 50. Scattering experiments jl4lL Il42j with 
isoscalar cr-projectiles have observed enhanced strength 
at E x « 0A2Eivgdr and its gamma-decay directly to 
the ground state indicating an isoscalar mode [19] which 
we label as low energy pigmy in Table I. For a globally 
applicable quantification the yield not seen as individ¬ 
ual spectral lines in a Ge-detector (looking like quasi¬ 
background) has to be taken into account 123]. 

For most of the nuclides to be discussed in the following 
( 146 Nd to 238 U) photon strength information for.E 7 < S n 
was obtained from individually known branching ratios 
of gamma transitions following neutron capture via res¬ 
onances near S n 53 or by analyzing primary ga mma 
spectra following average resonance capture (ARC) [125| 
to reach the nucleus in question. The absolute photon 
widths can be derived from 7 /n ratios and the reso¬ 
nance widths determined by neutron time of flight (n- 
tof). By an inspection of the gamma-ray angular distri¬ 
butions A = 1 is assured, if polarized neutrons are used 



FIG. 12. (Color online) Phot on s trength derived by photo¬ 
neutron production on nat Ba jllll ]f black circles) in compari¬ 
son to TLO for the IVGDR for 138 Ba (see the caption of Fig. [5] 
). Photo n sc attering data for E~, < S„ are sho wn as blue x- 
symbols Ill7tl and magenta stars for 136 Ba [l22ll reduced by a 
factor 0.7. 


the decay multipolarity (El or Ml) is known experimen¬ 
tally [95|. Inserting these widths and the average level 
spacings into Eq. right) results in /a( 15 7 ), but this 
relies on a known level density 1 /D r , and a ccor ding to 
Eq. (fTT!) this has lead to an increase by 30% [l43l ] in the 
respective table of RIPL-3, which is used in Fig. [3D] for 
the comparison to predictions on the basis of TLO. The 
dipole strength from such sources appears in the Figures 
as an addition to the IVGDR tail, suggesting that strong 
yield above it was not observable at other energies. 

In Fig. [13] the case of 146 Nd is shown and the regard 
of triaxiality in TLO leads to a reasonable description 
of data in the region of the IVGDR as well as below. 
In accord to Eq. CD and m comparatively small 
Ti of 2.82, 3.33 and 3.76 MeV are used without a de¬ 
crease with Ej. In previous work (3ll . mill a single 
Lorentzian (SLO, k = 1) was proposed for 146 Nd together 
with T ivgdr = 5.74 MeV, also shown in Fig. [l3]as red 
curve. It indicates that such a fit leads to a large reso¬ 
nance width T, as it emphasizes the pole region. Only 
with a decrease of the I VGD R width with photon energy, 
as was assumed [3ll [95l 1 144 ll45l | in KMF-type analyses 
[3l| . agreement to the low energy data near 5 MeV can 
be reached in spite of the large value for T. At even 
lower energy a component is added to the Lorentzian for 
the 7 -decay calculation which violates the Axel-Brink hy¬ 
pothesis. In 144 Nd and 146 Nd the KMF-model was shown 
[95l ] to produce a similarly good agreement to data like 
TLO as proposed here ( cf. Fig. [13]), but the latter does 
not need any modification of the Lorentzian shape. And 
it works well also for deformed nu clei , what is not the 
case for KMF in 163 Dy, as found [l44j ] from an analy¬ 
sis of 2-step gamma-casc ades afte r neutron capture. The 
data from 143 Nd(n, 70 ;) fl46l. [l47ll with their continuous 
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FIG. 13. (Color online) The photon strength for 146 Nd in 
comparison to TLO for the IVGDR. (dashed line in black). 
The data above S n are from photo-neutron production [lid ] 
and the ones below (magenta diamonds) are derived from 
gamma decay subsequent to ARC by 145 Nd [12511 . The dotted 
line in magenta results from a fit to these data with one pole 
only [za and the comparison to TLO is depicted as blue line 
like in the caption of Fig. [8] 

strength observed [l46 ] below 1 MeV are often quoted 
as experimental support of the KMF theory - negle ctin g 
the fact, that there Ml radiation may be f avor ed |l48| . 
similar to what was indicated for 94-96 Mo |l32j . KMF- 
theory was modified to improve the fit to photon strength 
data for various nuclei [125| . many of which are shown 
here to be in accord to our global approach. 

The data for 156 Gd shown in Fig. fill were selected as 
an example for experimental uncertainties to be aware of 
in discussions about details of the dipole strength: Two 
absorption experiments using rather similar techniques 
[l49l Il50j | disagree in the central IVGDR region beyond 
the ‘error bars’ published - not only in absolute but also 
in relative height for the two peaks, what does indicate 
an origin away from normalization or atomic background. 
No third data set is available and TLO comes closer to 
an average of both publications, when a beam energy 
spread of 0.3 MeV (rms) is applied. The apparent dis¬ 
agreement of the two datasets covering the IVGDR peak 
can be viewed as an indication of the reliability of the 
(q,n) data, many of which were obtained nearly 50 years 
ago. The ARC data included are subject to large Porter- 
Thomas fluctuations, which have been partly averaged 
out by re-binning the plots. The acceptable agreement 
found when comparing the data to TLO, as can be seen 
in Figs. [Q] [15] and H is in contrast to the strong ex¬ 
cess of the Lorentzian tail above the low energy dat a 
as reported previously for ‘non-deformed’ nuclei jl25| . 
Again it is stressed here that TLO uses smaller (photon- 
energy independent) damping widths leading to a smaller 
strength in the tail region. As discussed in previous work 
[13[ ]72{ a width dependence on excitation energy brings 



FIG. 14. (Color online) Photon strength for 156 Gd in com¬ 
parison to TLO (see the caption of Fig. [8]). Data are derived 
from ARC for 6-8 MeV (magenta diam onds) and from the 
cross section of photon absorption f [l49| ] black x-symbols and 
]l50ll black circles), which do not fully agree in the peak re¬ 
gion. 


the tail down considerably, causing a significant effect es¬ 
pecially in the region below 9 MeV (see Figs. [TT)landl24l 
and Figs. 1 and 3 of [H[). Within TLO this is avoided 
by having T depend on E » only and Eq. (fl2l) causes the 
two upper poles (see figs. [3 to 0, 0 andl2ll (and their 
sum) to have reduced height albeit all components have 
equal strength. Full conformance with the TRK-sum - 
is indicated near closed shells as well as for nuclei with 
large |Qo|- In Fig. [I5]results of two different experiments 
on 168 Er for the IVGDR range are depicted: The data 
with the larger error bars were obtained [H3 by photon 
absorption with subsequent subtraction of the strongly 
dominating absorption by the atomic shell, which has to 
be determined by a precise calculation. As the other data 
[l5lj stem from a photo-neutron experiment performed 
with a natural target, deviations might not be related to 
168 Er. But still the agreement of both to TLO confirms 
our parameterization, which does not yet include ISS dis¬ 
cussed with Fig. [G] The photon strength as observed in 
photon scattering (figs. [3 to [3 [3 . 13 and 0) shows 
intermediate “pigmy” structure at E 1 ~ 0.43 and/or 
0.55 Eivgdr, if these energies were covered experimen¬ 
tally. Also the q-decay observed after n-cagture is rea¬ 
sonably well accounted for in TLO (figs. 0 to 0), if 
an enhancement in the pigmy energy range is accounted 
for as minor strength, as parameterized in Table |T] Its 
contribution to the summed strength in Eq. m is 
weaker by at least one order of magnitude as compared to 
the IVGDR sum. Already many years ago, experiments 
with tagged photons [l 15t ] have identified a resonance¬ 
like structure in the cross section of photon scattering 
on targets in the vicinity of 208 Pb, and the case of Hg is 
especially significant. In 208 Pb the energy of such minor 
modes lies in a region of small level density and hence 
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FIG. 15. (Color online) Photon strength for 168 Er derived 
from ARC data (magenta * symbols) and from the cross sec¬ 
tion of photon absorption (black dots [l5Clt ] 1 in comparison to 
TLO (dashed curve in black). Shown as full line in blue is 
the prediction with minor strength and photo-neutron data 
for nat Er( 7 ,xn) [l5lll are depicted as black x-symbols. 



FIG. 16. (Color online) Photon str engt h for 190 Os derived 
from ARC data (magenta diamonds, [1521 ] 1 and fro m th e cross 
section of photo-neutron production (black dots Il53ll l. For 
comparison the TLO result is shown as before (see the caption 
of Fig. El). 



FIG. 17. (Color online) Photon strength for 196 Pt derived 
from ARC data Il25ll (blue diamo nds) and from the cross sec¬ 
tion of photoneutron production I154H (black dots) in compar¬ 
ison to T LO (see the caption of Fig. [8]). Photon scattering 
data [l22ll (magenta x-symbols) are re duce d by 30% in view of 
a ’’new” ansatz for the level density [124H . 



FIG. 18. (Color online) Photon strength in 202 Hg calculated 
by TLO for the IVGDR (see the caption of Fig. [H] ) is com¬ 
pared to photon scattering (magenta di amon ds) 115}] and 
photoneutron production (black circles) I155H . respectively, 
both studied with nat Hg. 


very large spacing between 1“-levels and it is intrigu¬ 
ing to rela te t his structure to the strong l _ -level at 5.51 
MeV [H, Il56| and eventually weaker ones seen nearby 
in this region of low level density p{E x ). This rather 
small p{E x ) reaching up to the IVGDR range is related 
to the large shell correction in near-magic nuclei as will 
be discussed in section ns Thus Porter-Thomas fluctu¬ 
ations have a strong influence on the data not leading 
to smooth strength function curves as assumed for the 
approximate pre dicti on us ed h ere. For 208 Pb several in¬ 
vestigations [83, Il55l Il57l ll58j have explicitly looked at 


fine-structure at E x < 12 MeV and some fluctuations in 
excess of the smooth IVGDR slope are seen; these are 
partly smeared out in the figures. Indications of contri¬ 
butions from the GQR are also ob served; unfortunately 
related electron scattering data do not allow 

a fully con sisten t transfer of information. Inelastic proton 
scattering il6l| indicates a peak at 21.5 MeV which could 
be either IVGQR or IVGMR, whereas the ISGQR is indi¬ 
cated in figs. IT 3 I to ll~9l. to partly overlap one or two of the 
IVGDR components. This excess over TLO contributes 
less than 5% to Eq. (ITOl) similar to the pigmy addition 
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contributing not more. Corrections of photon scattering 
yields have been performed using inelastic photon scat¬ 
tering data taken at the Urbana microtron laboratory 
[IH, where also the neutron emission directly to low 
levels of 207 Pb was observed [l58l |: it is of minor impor¬ 
tance for the arguments discussed here. In Figs. EH] and 



FIG. 19. (Color online) Photon strength for 208 Pb derived 
from ph oton scattering data using a quasi-monochromatic 
beam Ill5lb blue diamonds), from bremsstrahlung 
0 (magenta star symbols) and from two cr oss s ection 
measure ment s for photoneutron production [l55ll (black 
circles), fl57T (black x). TLO for the IVGDR is depicted as 
described for Fig. [ 8 ] 

|2T1 three sets of experimental data for 232 Th and 238 U, 
respectively, in the range of the IVGDR are displayed 
together: The data with the large error bars were ob¬ 
tained by photon absorption with subsequent subtraction 
of the dominating atomic absorption, which had to be 
determined by a precise calculation. They agree within 
uncertainty to data stemming from a photo-neutron ex¬ 
periment performed at Saclay, which were reduced here 
by 10%, as explained before. The agreement is not per¬ 
fect, but indicates the reliability of both in the IVGDR 
regime. The agreement between two data sets is impor¬ 
tant with respect to the disagreeing data obtained at Liv¬ 
ermore Il62j . These cross sections for 232 Th and 238 U 
are exceptional large in the sense, that an analysis on 
the basis of Eq. © indicates an overshoot of ~ 30% 
as compared to the TRK sum. One or even two pigmy 
modes seem to be present in photon scattering by all nu¬ 
clei with A > 70 regarded here. For lower energies the 
scissors Ml mode was predicted E3 to also be a general 
feature in all heavy nuclei, but in photon scattering it 
is difficult to separate from the El strength originating 
from quadrupole-octupole coupling. Results for a st rong 
presumably magnetic mode presented recently Il65| suf¬ 
fer from the missing parity assignment of the strength 
observed as well as from a questionable choice for the 
IVGDR tail. An identification as Ml is well possible by 
comparing e- and 7 -scattering mi and respective results 



FIG. 20. (Color online) Photon strength for 232 Th derived 
from photoneutron production cross se ction s from Saclay 
[l5i (black circles) and from Livermore [1621 ] fblack x sym¬ 
bols); both include photofission. Photon absorption data 
[16311 fblack + symbols) disagree to the latter, but agree to 
TLO (see the caption of Fig. [ 8 ] ) 



FIG. 21. (Color online) Photon strength for 238 U predicted 
by from TLO for the IVGDR (see the caption of Fig. [ 8 ] ) 
in comparison to two differing sets of data deriv ed from the 
c ross section of photo-neutron production Il55ll (black dots 
), [l62 ] (black x symbols); both include photofission. Also 
shown are photon absorption results [l63j ] (black + sym¬ 
bols with large error bars) and data from photon scatter¬ 
ing mi| (magenta ^-symbols) using a quasi-monoch roma tic 
beam. Data from 7 -decay after deuteron-scattering ll65ll at 
energies below 4 MeV are depicted as blue diamonds. 


indicate some difference to da ta fr om gamma decay after 
inelastic deuteron scattering Il65| suggesting a need for 
further study. For all nuclei presented here a quite good 
agreement between observations and the TLO prediction 
for El (plus minor) strength is found. TLO relies on the 
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classical sum rule (TRK), the IVGDR mode energies Ei 
are derived from the droplet model and their width Tis 
varying with it and hence slowly with A and Z. There is 
no dependence of it on E 1 , but on the three pole ener¬ 
gies Ei, which are derived from the deformation depen¬ 
dent axis parameters determined by CHFB calculations 
ll|. There is one deficiency apparent in these: In nuclei 
with or near closed shells the calculations predict more 
deformation 0 as deduced from B(E2)-values and, as 
found within the study presented here, in disaccord to the 
deformation induced splitting of the IVGDR. For closed 
shell nuclei a small deformation (reduced by a factor 0.4 
to 1 as discussed together with Fig. © results in a better 
description of the IVGDR peak shape; this reduction has 
no significant influence on the strength in the tail region, 
which will be shown to be essential for radiative capture. 
It is worth repeating here that the experimental data be¬ 
low Sn for isotopes used as a n ar gument for parameter- 
izations other than TLO [H, Il67| can well be described 
by TLO, if the extra ’minor’ strength is added. Also an 
account for the variance of the deformation parameters 
by instantaneous shape sampling [90 ] only leads to small 
changes of the calculated strength in the IVGDR for de¬ 
formed nuclei: The resonances are widened near the peak 
region, but the low energy tail remains unchanged. In 
closed shell nuclei and their neighbors the low level den¬ 
sity induces narrow variations in the cross section data 
up to energies near S n . These are e ventually identified 
as one type of pigmy resonances |168j; their contribution 
to the summed cross section stays below a few per cent. 
In general, the following conclusions can be drawn from 
the comparison of experimental data for 20 nuclei in the 
mass range from 50 to 250 to Eqs. CD and m- 


(a) The centroid IVGDR ene rgies derived from droplet 
model fits to masses 0 Il69( are in accord to the 
data, when an effective mass m e //c 2 = 800 MeV and 
the neutron radii as predicted from the CHFB calcu¬ 
lations [H are used. 

(b) There is no indication of a strong departure from the 
classical dipole sum rule 0, (54[ if the yield at en¬ 
ergies above the IVGDR is assigned to the quasi- 
deuteron strength. Previously reported excess may 
be due to overestimated resonance widths resulting 
from an ad hoc limitation to spherical or axial sym¬ 
metry. 

(c) The resonance widths vary only smoothly with A and 
Z and depend on the resonance energies via a power 
law (c/. Eq. (fl2l) h predicted by hydro-dynamical 
considerations [74|, I85 h 87I ] and supported by experi¬ 
ment Pd ]. 

(d) The IVGDR data together with data obtained with 
different methods for lower energies do not allow for 
a strong decrease of the width with photon energy, 
as previously postulated [Hj. Higher width values 
0 m 0 compensate such a decrease and cause 
large strength at low energy. 


(e) Only by allowing broken axial symmetry, three rather 
narrow resonance parts add up to a wide structure 
with a satisfactory large breadth of the IVGDR; this 
indicates the importance of calculations jlj or other 
input 0. 

(f) Experimental data may suffer from missing exact in¬ 
formation on detection efficiency as well as on beam 
resolution, eventually causing discrepancies between 
measurements at different laboratories- as was ob¬ 
served for 232 Th and 238 U . 

(g) The full recognition of the Axel-Brink hypothesis and 
the TRK sum rule together with the independence of 
the IVGDR width from the photon energy are fun¬ 
damental for our TLO-ansatz. 

The comparisons of a very similar approach to data for 
a number of is otop es as published previously 0, [0 HH, 
[H, [H, H3, HH, Ill4j ] can be considered part of the present 
work, such that a quasi complete sample of IVGDR data 
in heavy nuclei are shown to be well described bjyTLO. 
A description of the IVGDR in deformed nuclei |l70l ] not 
using the TRK sum rule and not based on a fully self- 
consistent calculation of the shape parameters was by far 
less successful in its predictions. The global TLO ansatz 
differs from locally adjusted Lorentzians presented pre¬ 
viously 0, [H, ll06( | and thus leads to other conclusions 
concerning eventual excess of strength in the tail of the 
IVGDR. 


VII. SPECIAL CONSIDERATIONS FOR ODD 
NUCLEI 

Concerning neutron s -wav e capture on J = 0 tar¬ 
gets it has been shown [l24j that the TLO ansatz also 
works for the odd nucleus 89 Y. When an average over 
the axis lengths for even neighbor nuclei is used to ob¬ 
tain / ei , a good description for the IVGDR region is 
obtained. It was assumed, that for Jo ^ 0 photon ab¬ 
sorption into a mode A populates m members of a mul- 
tiplet with m=min(2A+l, 2Jo+l) and the decay widths 
to the ground state Tor are equal for each member of the 
multiplet; the conditions for the validity of Eqs. (fldl and 
H are thus fulfilled. The strength observed corresponds 
to the cross section summed over the multiplet and this 
can be described by an effective g , which according to 
Eq. (© is: 


9 eff = 


E 


2 J r -f-1 
2J 0 + 1 


= 2A + 1 


(15) 


This ansatz is valid in heavy nuclei 0] as it relates to 
the condition of weak coupling between the odd particle 
and the mode A. The TLO-calculations for odd-A nuclei 
as shown in figs. 0 to 0 were performed on the basis of 
Eas. (fTTirH21 and[T3l) with k = 3. No extra spin dependent 
factors are needed and agreement to the experimental 
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FIG. 22. (Color online)Photon strength for 133 Cs calculated 
as described in the text (see the caption of Fig. [8] ), in com¬ 
parison to tw o da tasets of photo-neutron production shown as 
black points [l37j | (like all data from Saclay reduced by 10%) 
and black x-symbols Il5ll . The photo-absorption data below 
Sn (magenta asterisk) are de rive d from elastic scattering by 
the neighbour nucleus 132 Xe Il23ll . modified by 0.7. 


data is found to be similar as in the case for even nuclei, 
also in the tail region below S n . In Fig.[55]data for 133 Cs 
are shown to be close to those for neighbouring even nu¬ 
clei depicted in Figs. fTSl and fl3l The agreement to TLO 
is obvious and also “minor” strength as pigmy and the 
IVGQR are seen. For the nucleus 197 Au not only the 



FIG. 23. (Color online) Photon strength derived for 181 Ta 
from photon scattering ll23H fmagenta ^-symbols, below 7 
MeV, modified by 0.7) and from photoneutron production 
(black x symbols, [l 71 1 ] : black dots, f 1~5 it ] ) in comparison to 
TLO (see the caption of Fig. [8]). 

prediction is depicted, but also a SLO curve from RIPL- 
3 [HUzzI, which clearly over-predicts the data below S n 
extracted from Fig. 18 of ref. (l5j| . In contrast to the 


missing strength as compared to a single Lorentzian used 
there, the agreement is reasonable for TLO, when the 
discontinuity near 19 MeV is related to the known [108} 
incorrect separation of the 2n-channel. The inclusion of 



FIG. 24. (Color online) Photon strength for 197 Au derived 
from photon scattering data ' 15j (magenta asterisks) using a 
quasi-monochromatic beam and from the cross section of pho¬ 
toneutron production [l55ll (black circles) in comparison to the 
TLO prediction (see the caption of Fig. [8]). Also shown are 
newer data for the peak region 1 1041 (black x- symbols) and a 
SLO fit curve from RIPL-3 (magenta dotted line). 

triaxiality in TLO leads to a reduction of and thus of 
a a bs for sufficiently large (E? — E 2 ) 2 in Eq. (fill) . The 
widths Ti used previously JT^j are 2.9 and 4.0 MeV and 
an additional factor of 1.22 was obtained as compared to 
the TRK-sum rule. This factor is 1.0 for TLO and the 
values for are 2.7, 3.0 and 3.5 MeV. When 197 > 198 Au 
is considered spherical T « 4.5 MeV results from a SLO- 
fit and this eff ects largely the strength predicted in the 
tail region [H, 150}. Also for 197 Au the use of the KMF- 
model was proposed [95} and the more satisfying agree¬ 
ment of TLO as presented in Fig. [24] again favours our 
ansatz over this model, as already observed for 146 Nd, 
Fig- EH In the bombardment with photons 239 Pu mainly 
undergoes fission and the weaker neutron emission chan¬ 
nel has to be added to obtain a a b s - In Fig. [55] the result 
of this sum of Livermore data jl72j is compared to direct 
absorption data and a reasonable agreement is seen, as 
well as a good agreement to the TLO-prediction. This 
is remarkable in view of the disagreements depicted in 
Figs. M and [5l] for the near neighbors 232 Th and 238 U 
and doubts about these older data seem justified. 

Together with the agreement between 88 Sr and 89 Y re¬ 
ported recently |l24j the examples presented in this sec¬ 
tion support the TLO ansatz for the derivation of photon 
strength in odd nuclei. Our assumption, that the damp¬ 
ing widths T; do only depend on £) and not on E 1 clearly 
is at variance to previous_proposals made for the electric 
dipole strength function [95} and recently made [3lj| mod¬ 
ifications of this work. But we will show, that our new 
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FIG. 25. (Color online) Photon strength for 239 Pu (black 
dashed curve: TLO; full blue curve: minor strength added to 
TLO). Data are derived from summing cross sections of fis¬ 
sion and neu tron emission induced by quasi-mono-energetic 
photons Il72l1 (black circles) and by discrete 7 -rays from neu¬ 
tron capture [1731 1 (magent a st ars on the low energy slope). 
The absorption data from Il63il (black x) were obtained with 
bremsstrahlung. 


ansatz allows to calculate average radiative strengths and 
neutron capture cross sections which compare well to ob¬ 
servations. 


VIII. PHOTON STRENGTH OF OTHER 
CHARACTER THAN ISOVECTOR ELECTRIC 
DIPOLE 


The height of the low energy tail is nearly proportional 
to the IVGDR width and it depends weakly on its de¬ 
formation induced splitting. A consequence hereof is a 
nearly full independence of the strength on any free pa¬ 
rameter. But literature on photon scattering in this en¬ 
ergy range has pointed out the possible importance of 
strength of other character than isovector electric dipole, 
as recently reviewed M- In the following also work pub¬ 
lished more recently is included and the discussion will 
be generalized to electromagnetic strength of different 
character. Even if it has minor importance for photon 
absorption and the sum rules, it may influence the decay 
of the compound nucleus and consequently also the re¬ 
action cross sections. Strength eventually adding to the 
IVGDR tail is characterized by approximate expressions 
listed in Table U and its eventual importance will be in¬ 
vestigated. A direct relation exists between ground state 
transition widths, summed within an interval A e, and 
the /ex (in GeV“( 2A+1 )) as derived from Eos. (THU fill) . 
Strength information can hence be obtained by summing 
spectroscopic width data (in MeV) over a given energy 


range A e (also in MeV): 


1 




■ r 7 


a b ^ El x+1 


T x 

abs 


(:irhc) 2 g e ffE 2A X A E 


(16) 


The second part of Eq. (fl6l) allows to relate the strength 
function to integrated absorption cross sections I x bs ; 
these can be directly related to sum rules. 


Electric dipole strength below the IVGDR 


A low energy dipole mode was predicted to be formed 
by E3 strength coupled to low energy quadrupole modes 
11]. In many even-even nuclei rather strong photon ab¬ 
sorption into 1 “-levels with IMeV < E x < 4MeV has 
been observed 1174 ] and a correlation of E x to the sum 
of the excitation energies of the low collective 2 + and 3“- 
modes was established. A systematic dependence of this 
dipole strength on A and Z has not been registered, prob¬ 
ably due to possibly destr uctive interference between dif¬ 
ferent contributions [l75| . A functional dependence pro¬ 
posed for the 2 + x 3“-mode from theoretical arguments 
ta was further specified GMz! and the B(El)-values 
used here are in the range of observations in nuclei of 
various degrees of deformation. Photon scattering stud¬ 
ies show similar strength for odd and even nuclei [174| . 
but fragmentation away from closed shells and especially 
in odd nuclei cannot be excluded; it may cause a missing 
of small strength components. To estimate the influence 
of this low energy dipole mode on the photon emission 
we approximate it by a Gaussian distribution with pa¬ 
rameters listed in Table [T) The centroid energy E qo is 
taken as the sum of the energy of the 1 st 2 + -level [I 4 ] 
and th e cor responding value for the octupole (3“) exci¬ 
tation fl77| : the oret ical approximations for exotic nuclei 
are available QLUzi- 

In principle Eq. m allows the inclusion of numer¬ 
ous recently published results for strong peaks not far 
below S n near 7 MeV. In case an assignment to El is 
made, this strength is often denoted as PDR (‘Pygmy 
Dipole Resonance’) [l9j, Il79| . apparently in analogy to 
Bartholomew’s “pigmy resonance” in strength functions 
pbj . Such bumps were assigned to the “vibration of ex¬ 
cess neutrons against a proton-neutron core” [l80j . or 
also as “an integral part of a toroidal El mode represent¬ 
ing an example of vortex collective motion in nuclei” [l79( | . 
“Intermediate structure” in photon scattering was dis¬ 
cussed in detail [80] and the concept of photon strength 
functions f\{E 1 ) was then introduced. At energies ap¬ 
proaching S n from below an extraction of f\(Ej) from 
photon scattering data becomes questionable as compet¬ 
ing channels suppress the emission of elastic photons and 
often unobserved decay branches (inelastic photon scat¬ 
tering) are not easily accounted for. This may fake a fall- 
off which, together with the rise when approaching the 
IVGDR, appears as a peaking distribution. The strength 
in that energy range identified in the isotopes 112_124 Sn 
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El El by high resolution 7 -spectroscopy was used 
to estimate the full dipole strength in the “PDR” by 
accounting for additional quasi-continuum contributions 
and the loss due to branching. Its peak energy increases 
weakly with decreasing A, but a clear change of the in¬ 
tegrated strength is not observed even when combining 
the observations with data for 138 Ba, 140 Ce |l42l ] and 
recent work for Xe-isotopes El- The analysis of ob¬ 
servations on seemingly isolated peaks in gamma spec¬ 
tra deserves extra care: Porter-Thomas fluctuations may 
sham structures in data taken with a resolution larger 
than the level distance and thus not sufficient to resolve 
all levels jl 8 l| : Sampling will give many small numbers 
and a few large ones and this is ideally suited for pro¬ 
viding the pronounced peak structure that makes for the 
happiness of every dyed-in-the-wool spectroscopist. 


In a few earlier papers this problem was addressed ex¬ 
perimentally with special care in the photon detection 
and strength data were published for structure seen in 
the region between 5 and 6 MeV with photon scatter¬ 
ing by Zr and Sn [l| as well as for A ss 200 |l!5j . As 
mentioned there, the contribution of the quasi-continuum 
below the lines is significant even after the real back¬ 
ground due to unwanted radiative processes in the de¬ 
tector and the near-by environment has been quantified 
and subtracted. In the limit of E 1 —> 0 the strength fEi 
shold vanish as suggested by the well-known fact, that 
nuclei do not have an electric dipole moment S®. Un¬ 
fortunately Thomson scattering dominates at long wave¬ 
lengths and obscures resonant scattering processes. Frag¬ 
mented parts of the electric dipole strength outside of the 
IVGDR may have isovector or isoscalar character and a 
distinction by exper iments with isoscalar beams was pro¬ 
posed [191 ll4ll i142( | . This approach is not unambiguous, 
as only in collective models based on one incompressible 
fluid isoscalar electric dipole strength is isospin forbidden 
as it would cor respo nd to a spurious motion of the cen¬ 
tre of mass 11, E|- Hints for the isoscalar character of 
electric dipole strength may indicate non-uniform proton- 
neutron distributions or compressional modes. Strength 
between E x = 5.5 MeV and the neutron separation en¬ 
ergy S n was shown to be of isoscalar nature in 40 Ca, 
58 Ni, 90 Zr and 208 Pb by the coincident observation of in- 
elastically scattered a-particles and de-excitation 7 -rays 

El. 


In the present study the various results reported for 
this energy range are considered an addition to the TLO- 
extrapolation and tentatively separated into two compo¬ 
nents named low and high energy pigmy mode (PM) in 
Table IU the parameters given there are selected such that 
the total strength included may be understood as an up¬ 
per limit. In figs. 0 to El the low energy PM (a seemingly 
resonant strength near 0.43 x Eivgdr) appears to be of 
similar magnitude for many different A, when compared 
to the IVGDR. Hence this intermediate structure resem¬ 
bles a prediction made for 208 Pb: In a three-fluid hydro- 
dynamical model of heavy nuclei El an oscillation of 
the excess neutrons against a core with N = Z was shown 


to result in a peak of energy E x < 0.45Eo, with E 0 valid 
for the full vibration of the neutrons against the protons, 
the IVGDR [68l - l7| . Experimental “evidence for a 5.5- 
MeV radiation bump”, an intermediate structure then 
named ‘pigmy’ El, suggests our nomenclature high en¬ 
ergy PM. Independent of an interpretation of these PM’s 
valid for all heavy nuclei a Gaussian seems justified for 
both and the parameters used for them are given in Ta¬ 
ble n These values were selected to clearly overpredict 
most available data after the TLO-integral over the en¬ 
ergy interval is subtracted El- They are purely phe¬ 
nomenological and independent of isospin and of models 
like the ones proposed ® to eventually explain extra 
strength below S n . 


Magnetic Dipole Strength 


For both dipole modes (El and Ml) Eq. (TTBli can be 
used together with the generalization of Eq. © to gather 
strength information from spectroscopic data. Here, the 
squared reduced matrix elements are inserted, which are 
direction independent in contrast to the still often used 
B(E, M A) values. They are, expressed in fm 2 , obtained 
from: 


r r0 (£ 7 ;El,Ml) — 


167T a e E 3 
9 g(hc ) 2 


(r||El, M1||0) | 2 (17) 


The photon strength between nuclear ground states and 
the IVGDR region is predominantly of electric dipole 
(El) character and also in the tail region around 5 MeV 
Ml decays are weaker as compared to El [95]. Magnetic 
(Ml) spin flip strength occurs at higher energy than col¬ 
lective orbital magnetic strength (scissors mode), which 
is strong in nuclei with large quadrupole moment and 
possibly contributes to radiative capture. Both are built 
either on ground states or on excited le vels a nd then con¬ 
tribute to the strength on top of them [ 1831 ] in the sense 
of the Axel-Brink hypothesis. High r esolu tion photon 
scattering data as e.g. published by El show spec¬ 
tral details indicating low level density, and these are 
especially significant in near magic nuclei. Again weak 
components at low energy are likely to be hidden in a 
quasi-continuum partly due to experimental background. 
Combined measurements of transition rates and polariza¬ 
tion in the gamma decay of compound nuclei have been 
performed e.g.for A = 79 jl84l | and the strength of low 
energy Ml transitions has been determined to be below 
that for El. But measurements are rare for statistically 
significant information on El and Ml transitions in one 
nucleus. A possible way out is to a verag e over many nu¬ 
clei and a suitable compilation [76L Il85| is available and 
updated regularly. The sorting of nearly 3000 transitions 
from compiled data ® for 80 < A < 180 has resulted 
in the plot presented in Fig. [26] The Ml-strengths 
/m i(-E 7 ) are apparently considerably larger than the El- 
strengths /ei(E 7 ) for E 1 below 1 MeV whereas above 
the El strength becomes stronger up to the region of the 










22 



FIG. 26. (Color online) Strength functions for an ensemble 
of nearly 3000 transitions in heavy nuclei with 80 < A < 180 
as extracted from a data compilation 51 EH- The vertical 
bars indicate 20% of the root-mean-square of the distribution 
of the data in each bin and the points depict the averages for 
Ml (blue diamonds) and El (black circles). 


IVGDR. A direct measurement of the transition strength 
contributing to the decay width in the high level density 
region near S n was possible by using experimental infor¬ 
mation gained in 143 Nd(n, 7 a) experiments mm - in 
view of Fig. [211 this “low energy” stre ngth is likely to be 
mainly Ml as pointed out before EH , but El was occa¬ 
sionally assumed to be the main component [95j . In any 
case the summed dipole strength function /i in this mode 
was shown 132L [l46j to be small and concentrated at very 
small E 1 . Although such a low energy Ml-component 
has been used [95), ll86l | eventually, we consider it unim¬ 
portant for the prediction of radiative n-capture: The 
compilation of experimental data shown in Fig. dH de¬ 
picts a significant decrease for E-y > IMeV and for lower 
energy the cross section will be shown to be suppressed by 
the level density in the final nucleus. Recently very low 
energy Ml transitions resulting from orbital rearrange¬ 
ment have be en p redicted by shell model calculations for 
Mo-isotopes [l32| to eventually compete to El below 3 
MeV. 

In a recent review 53 the properties of three compo¬ 
nents of magnetic dipole strength have been described 
in detail: For energies between 2 and 4 MeV orbital Ml 
strength has been observed, and it has been shown that 
with increasing ground state deformation f3 this orbital 
scissors mode reaches an integrated strength proportional 
to Z 2 P 2 at an excitation energy 0.21Ejvgdr, in agree¬ 
ment with data reviewed recently [1871 ]. The enhance¬ 
ment of gamma decay between high lying levels as ob¬ 


served in two step cascades [1281 . Il45j was attributed to 
orbital scissors mode of about the strength seen by pho¬ 
ton scattering. In addition to it a double humped struc¬ 
ture was seen at higher energy (around 41 A -1 / 3 MeV) 
and identified 53 as isoscalar and isovector spin-flip Ml 


modes. This is in agreement to the strength observed 
in heavy nuclei with polarized neutrons [951 ]. Including 
both regions at ~ 34A -1 / 3 MeV and ~ 42A -1 / 3 MeV 
53 the energy weighted sum can be calculated within 
the energy range A e with reasonable accuracy; it reaches 
< 2 MeV fm 2 , as the A-dependence of peak energy and 
peak integral approximately compensate each other 53- 
The resulting estimate is smaller by at least 2 orders of 
magnitude as compared to the main component of the 
GGT sum rule (eq. 5.6 of [H}). This prediction was 
derived from very general arguments for the total pho¬ 
ton strength in heavy nuclei, i.e. including all multipo¬ 
larities. In Table Q] the estimates used by us are given; 
the scissors mode is considered as well as the two spin- 
flip Ml components, which are becoming important at 
higher photon energy. Polarization data are necessary 
to identify magnetic strength unambiguously a nd o bser¬ 
vatio ns m ade with polarized photons on 88 Sr 513 and 
90 Zr [l88l ] were performed at two different laboratories 
and with two different techniques: off axis continuous 
bremsstrahlung and tagged photons. A surprisingly high 
magnetic st rengt h was recently used in calculations for 
Zr-isotopes |l05j . although it superseded the systemat- 
ics 53 by a factor of ~ 3. Combined to a rather low 
El-s trength as derived from a HFB+QRPA calculations 
[l89j [ a suitable agreement to capture data results; in view 
of the Ml-enhancement needed the underlying El model 
may have to be questioned. 


Electric quadrupole modes 


Low energy E2- transitions can be considerably en¬ 
hanced and have played an important role in the spec¬ 
troscopy of heavy nuclei - as discussed in sections H3 
andlllJ. The photon absorption cross sections can be de¬ 
rived from Eq. ©; they have been found to be small 
as compared to El absorption in the energy re¬ 

gion near S n . In figs. @ to 53 Je 2 was multiplied with 
E 2 , the additional energy dependent phase space factor 
for absorption and decay, to allow a visual comparison 
to /i; but for the calculation of the integral Ie 2 = Ic 
(as listed in Table Q]) the correct d ecay width i s use d in 
Eq. ([HI). Extending previous work [bl l3ll. fill Il90] the 
influence of quadrupole giant resonance (QDR) contribu¬ 
tions to photon absorption was investigated; information 
for quadrupole strength com e s fro m sum rules and the¬ 
oretical predictions Jl0l. Jll6l . EH adjusted to electron 
scattering data [H, Ho, Il92j . The isoscalar ISGQR lies 
not far from the pole of the IVGDR such that it adds 
to the strength observed there. Strength modifications 
due to the IVGQR appear to be important, as seen in 
figs. 53 to 53 , and result in a difference from TLO at 
higher energy. The influence of quadrupole modes at 
even lower energy is investigated by assuming 5% of the 
ISGQR integral J a abs dE to be downshifted by a factor 
of three in energy. The approximations used in Table U 
for the integrals of the El, E2 and Ml components were 
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Component 
parameter (units) 

multi¬ 
polarity E c (MeV) 

I c (fm 2 MeV) 

(MeV) 

low E x pigmy mode 

El 

0.43£o 

rj Z(N — Z) 

‘ A 

0.6 

high E x pigmy mode 

El 

0.55£o 

13 Z(N ~ Z) 

0.5 

0+ <H- (2+ x 3“)i_ 

El 

W ( 1 + ¥) 

0.006ZA/3 

0.6 

orbital (scissors) mode 

Ml 

0.21E 0 

0.033 ZA/3 

0.4 

isoscalar spin-flip 

Ml 

A2A~ 1/3 

17 

0.8 

isovector spin-flip 

Ml 

47A _1/3 

27 

1.3 

low E x quadrupole 

E2 

19A _1/ 3 

n 1 OL e (irR c Ery) 2 Z 2 
(3 Am p c z ) 

1.0 

ISGQR 

E2 

63A _1/ 3 

0 ct e (irR c Ery) 2 Z 2 
(3 ArripC 2 ) 

1.0 

IVGQR 

E2 

48A“ 1/6 

ct e (TrR c E y ) 2 Z 2 
(3 Am p c z ) 

1.8 


TABLE I. Parameters of an upper limit for three minor electric and three magnetic dipole modes to the dipole strength 
function to be calculated with a Gaussian fEa 1 1911 . A very similar shape near the peak is reached for a Gaussian with a rms 
breadth a, which is larger by a factor 2.5 as compared to the width P of a Lorentzian. Three GQR modes are also listed, which 
are seen in photon absorption as well; they are proportional to the square of the nucleus’ charge radius R c (rms). 


derived through a comparison to respective observations 
fl7l l95l. fl88l| under respect of the above mentioned in¬ 
formation. The resulting parameters for these “minor” 
contributions to the photon strength function are listed 
there; some of the components (e.g.Ml and ISGQR) have 
been enlarged to assure that an upper limit is used in the 
subsequent estimates for the influence of “minor” modes 
on photon absorption and emission. Separate Gaussian 
distributions: 


h(E y ) 


1 _^ ex / -(£ 7 -£ c ) 2 \ 

(jthc) 2 g e ffE^ x ~ x y/2Tt(T c \ 2(J c ) 


( 18 ) 

for each of them are used to avoid unobservable strength 
in long tails. No arguments were found, why Lorentzians 
should be preferred; this differs from the IVGDR, where a 
Lorentzian with energy independent width is appropriate 
[3 El]- The I c in Eq. (ITSl) as well as in Table Hare sums 
over many individual levels within the Gaussian width 
tr c like in Eq. m for the IVGDR and at variance to Eq. 
0 , valid for a single level. 

Table U lists the cross section peak integrals I c as well 
as the centroid energies E c and the breadths a c (in MeV) 
as discussed above and used for the approximate consid¬ 
erations presented in sections Hxl t olxll. In figs. @ to [3 the 
corresponding results are shown as blue lines. The rela¬ 
tive importance of the components listed is not easy to vi¬ 
sualize for the many nuclei under study here, as there are 
considerable variations with A and Z . As will be shown 
in the next section the level density increases strongly 
when approaching S n and due to the convolution with 
the photon strength function the sensitivity of radiative 
capture peaks at E x ss 3 — 5 MeV [83] - far below S n . 
Clearly a significant contribution to this process is due to 
the electric isovector dipole mode described by TLO. The 
other important contributions to the photon strength are 
related to the pigmy resonance and to the orbital mag¬ 
netic “scissors” mode , at least in nuclei with suffi¬ 
ciently large deformation (/3). The spin flip Ml strength 
values listed in the Table are upper limits deduced from 


experimental information H3, HHI, Il 88 ( | ; this mode with 
its elevated energy has a reduced influence on radiative 
neutron capture. In contrast, the electric dipole strength 
due to the phonon coupling of the type (2+ x 3“)i_ is 
of more importance; but it is difficult to quantify it in a 
scheme which covers the full range of nuclides in ques¬ 
tion, and thus our predictions clearly have some uncer¬ 
tainty. It also is mentioned, that for transitions between 
two high lying levels El, Ml and E2 transitions may be 
important that are not accessible by studies with photon 
or electron beams. 


IX. LEVEL DENSITIES 
Intrinsic state density 

Nuclear level densities p(E x ,J v ) determine the final 
phase space for predictions of compound nuclear cross 
sections and decay rates. Within the RIPL project [3l| 
various attempts were reviewed aiming for a simultane¬ 
ous fit to data for s-wave neutron capture mea n spa cings 
and low energy discrete levels. It was shown [1931 ] that 
a parameterization based on Fermi gas formulae (modi¬ 
fied to finite systems) leads to unsatisfactory predictions, 
but no mention of broken axial symmetry as possible 
cause can be found. In calculations based on fundamen¬ 
tal principles, of interest especially if they need few free 
parameters only, a clear distinction has to be made be¬ 
tween the intrinsic quasiparticle state density ui qp (E x ) 
and the level density p(E x ,J 7V ) in the observer’s sys¬ 
tem. For the calculation of io qp on an absolute scale we 
use an analytical approach. It avoids as far as possi¬ 
ble parameter adjustments and strongly relies on statis¬ 
tical laws for a Fermi gas - a system of independent spin 
1 /2-particles with weak mean attraction. As shown long 
agopjMj it is characterized by a gap A (f) falling with 
rising temperature parameter t down to 0 at a ‘critical’ 
t pt = AoexpC/ir = 0.567 • Ao (with the Euler constant 
C=0.5772). This indicates a phase transition of probably 
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second order |l95l4l97 |. governed by pairing. In infinite 
nuclear matter (nm) the backshift energy E bs , which de¬ 
scribes a shift between the Fermi g as ze ro and the nu¬ 
clear matter ground state, is equal |l96j | to the pairing 


condensation energy E con , defined below in Eq. |23l 


An application of text book statistical physics 11951 to 
nuclear matter allows to evaluate the general features of 
this phase transition, and all effects ap peari ng in finite 
nuclei can be trea ted micro-canonically |196j . Foll owin g 
earlier work jl98l ] and a recent data analysis |l93l . Il97j , 
the energy dependence of the state density is assumed 
to be exponential in the pairing dominated phase below 
the phase transition point. This region is hence approxi¬ 
mated by the assumption of quasi constant temperature 
T c t, to be distinguished from the Fermi gas temperature 
t). It will be discussed below, how these phenomenolog¬ 
ical parameters and 0 J qp ( 0 ) are found for the nuclei of 
interest. 

If only quasiparticle excitations are considered the to¬ 
tal state density (in the intrinsic frame) ui qp (E x ) at exci¬ 
tation energy E x is approximated by: 

io qp (E x ) = u qp ( O)exp(-A) for E x < E pt (19) 

J-ct 


^qp(E x ) 


s/k exp(2^/a(E x - E bs )) 
\Ta'l\E x - F fcs )5/ 4 


for E x 5^ E p t 
( 20 ) 


At the phase transition energy - corresponding to t pt and 
given by E pt = a nm ■ t pt + E bs - a transition from a Fermi 
gas like behavior above to a pairing dominated regime 
below the transition point occurs. For the Fermi gas 
phase of infinite nuclear matter (nm) the “level density 
parameter” a nm is required to be inversely proportional 
to the Fermi energy £f 0 ■ It also determines the e nergy 
Ernn of t he pairing induced condensation [H, IlSfil Il9fil 
1991 l200j . For (E > E pt ) in hnite nuclei, the backshift 
Eb s stands for the energy between the Fermi gas zero 
and the nuclear ground state as in nuclear matter, but 
it has to be corrected for shell effects. In addition a 
surface term 6a has to be added to a nm to obtain a level 
density parameter a for a given A. Using a in Eq. (1211) to 
quantify it, a comparison to capture resonance spacings 
finds a to be as small as a = 0.06, when treated as a 
parameter, globally valid for all A under regard - actually 
the only one in the present fit to neutron capture data: 


Q"nm 


^- = ^a = cinm + Sa; Sa = aA 2/3 (21) 

4 £f 15 


J7 — 3 A 2 

J-^con 27T 2 u ' nrnL ~ 1 0 

E bs = E con — 6E(Z 1 A) (22) 

E p t — crtpt A E^ 

Various approximations are required jl24 ] ; they are listed 
here albeit not all of them are of significant influence for 
the conclusions made later: 


1. The pairing parameter A((E X = 0) is approximated 
by A 0 = 12 -A -1 / 2 , independent of angular momen¬ 
tum. 


2 . Ao is used for neutrons and protons and thus inde¬ 
pendent of neutron excess N — Z. 

3. Quasi-particle states are evenly spaced (on average) 
at the Fermi energy, not varying with N — Z. 


4. Fermi energy £p = 37 MeV and nuclear radius R = 
1.16 • A 1 / 3 are independent of N — Z\ they control 
a nrn for nuclear matter and the phase transition 
[l95l . Il96l ] at excitation energy E pt . 


5. A direct dependence of equilibrium deformation on 
excitation energy E x and angular momentum J is 
neglected beyond its implicit influence via SE(Z , A) 
in Eq. (l23l) . 


6 . A factor 1/4 accounts for i? w -invariance [201 1. 


7. The influence of shell effects is controlled by 
6E{Z 1 A), extracted by subtraction of the experi¬ 
mental mass from liquid drop values, following ref. 
10 2 


8 . At variance to previous work p. [IM [M [ 200 I 
|203|| the shell correction ter m is direc t ly ap plied 
to the backshift energy E bs [ 198i . 2041 . 12051] . No 
dependence on deformation is included. 


The expression given by Eq. (l20l) for uj qp in the Ferm i 
gas regime - initially derived neglecting pairing [ 6 j, ll98l ] 
is in infinite Fermionic systems a good approximation for 
the f ormalism derived with a thorough (micro-canonical) 
[l96i | inclusion of pairing, if E x is back-shifted by the 
condensation energy E con . In finite nuclei the back-shift 
Ebs as given in Eq. G51) is the sum of the pairing term 
E C on, approximately independent of A, and the effective 
shell correction 5E(Z,A) which includes the odd-even 
mass difference. The intrinsic (quasi-particle) state den¬ 
sity uj qp (E x ) for the Fermionic region as well as for E pt 
are given by Eqs. [20] to [23[ Below E pt an interpolation 
of u}q P (E x ) to the ground state is found by Eq. (fl9l) : the 
simple approach of a logarithmic interpolation in analogy 
with an exponential inc rease of jj nr , ( E x ) is in agreement 
to observations [HI, IM [M EMail- 
We fix E p t using a, t p t and E bs and the requirement of 
a continuous transition in uj qp (E x ) at E pt to determine 
T ct . In principle w 9 P (0) can be fixed locally by regarding 
known spectral data at low (E x , similar to what has been 
done previously P, Il22l Il90j . In a less stringent way 
we used as a first global approximation just above the 
ground state w gp (0) = 0.1/A o . The corresponding effect 
is demonstrated below with the predictions for average 
radiative widths. Eq. JUD shows, that a temperature 
parameter t = \/(E x — Ehs)/a for a Fermi-gas approx¬ 
imated [!, Il96l Il98j at the saddle point of the Laplace 
transform, differs from an apparent nuclear temperature 
T app = as well as from T ct , which turns out to be 
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smaller than T app by up to 35%. The resulting sudden 
change in the slope of uj qp (E x ) at E pt {cf. Fig. [5HJ) is 
acceptable for a second order phase transition. Close 
to magic nuclei the large (negative) shell correction re¬ 
sults in a then large energy E pt for such a break. It is 
remarkable that in weakly bound nuclei ( e.g. near the 
edge of stability) the Fermi gas phase starts already be¬ 
low E x = 3 MeV. It can be seen in Fig. [271 that for 
more than 100 of the 146 nuclei investigated here E pt is 
smaller than S n and thus the neutron capture resonances 
fall into the Fermi gas regime, but the subsequent gamma 
decay mostly ends below E pt and the level density there 
is important for the prediction of radiative capture cross 
sections. In the Figure it also becomes obvious, that in 
the Fermi gas regime SE(Z , A) is closely correlated to E pt 
and thus also to uj qp {E x )-, it is hence a quantity of large 
importance for a level density prediction and S n plays a 
minor role, at variance to a previous assumption [198J. 



A 


FIG. 27. (Color online) Phase transition energy E p t for nuclei 
in the valley of stability vs. A (full line in black) in comparison 
to values for S n (dashed in magenta) and the shell correction 
energy 5E (blue dotted line), all in MeV. 

The quantities to be compared to observed level spac- 
ings have to be derived from ui qp {E x ) by a projection 
on angular momentum J in the o bser ver system. The 
proposal was made d, EM ES Hq 3 to consider the 
m-substate distribution of uj qp {E x ) as a Gaussian with 
width o' around m = 0 and to differentiate atm = J+l/2 
with respect t o to. This leads to a sp in dependent level 
density ffil. l8ll. H67L flSfiL H9fil. : 

fr , 77n „ 2J+1 ( (J+1/2H /c , ^ 

p sp h{E x ,J ) - 2v / 87T(7 3 exP 2 ct 2 J w 9p(^) 

(23) 

The spin dispersion a can be related to the nucleus’ 
moment of inertia 3, often assumed to be the rigid ro¬ 
tor value [|]. This redistribution of the quasi-particle 


smallj 2 J T 1 


2-\/87ro 


■w qp (E x ) with (7 = 


states into levels of distinct spin implicitly assumes [207] 
the nucleus to be exactly spherical symmetric even at 
E x — S n - This neglects strongly mixed modes which, 
due to their collectivity, are pulled from their original 
quasi-particle energy down into the low excitation region. 
Without assuring spherical symmetry for the nuclei stud- 
has found a widespread use [HI, mm 
|204[. In some work the rotational collectivity 
present in an axially symmetric nuc leus was included at 
this stage as extra term [a i20il20fiU207| . yielding a level 
density enlargement by a factor a 2 ( « 21/5) as com¬ 
pared to Eq. (17^1) . But still an agreement with observa¬ 
tions was not reached without a sign i fican t enlargement 
of a as compared to a nm [6j, HI EH, l206| ; often it was 
adjusted in a fit. 


Collective enhancement 


It is obvious H9i l20H that the absolute scale of any 
prediction for the level density is strongly influenced by 
breaking spherical symmetry in the Fermi gas regime and 
at E pt . In previous work [3l| deficits in the comparison to 
resonance spacings observed at the neutron binding en¬ 
ergy S n have been compensated not only by an increase 
of the level density parameter a [311, but also by an ex¬ 
citation energy dependence of it. Both measures are not 
supported by a ’’standard theory” of a Fermi gas, which 
we accept as proper description of the statistics in highly 
excited nuclei. We will show by a comparison to experi¬ 
mental data that major modifications of a are not needed 
when giving up the scheme valid for the case of spher¬ 
ical or axial symmetry. An inclusion of the collective 
enhancement of p by breaking axial symmetry follows a 
proposal made long ago [20 lj ] and also presented in text 
book of Bohr and Mottelson [f|, cf. Eq. (4-65b). One 
obtains from this an astonishingly simple expression for 
small J, which to our knowledge was not taken up in 
any subsequent comparison to data; a factor 1/8 for 
invariance and parity conservation leads to: 


p{E x ,J*)^™-±± Ugp (E x ) 


(24) 


In view of the triple split observed in the IVGDR of 
nearly all nuclei with A > 60 - as described above - 
it is quite obvious to account for triaxiality in the en¬ 
hancement of level densities. It results from the build-up 
of a rotational band on each intrinsic quasi-particle state: 
The total level spectrum, for a given angular momentum, 
is therefore obtained by summing over a set of intrinsic 
states rather than by a decomposition of the level spec¬ 
trum, as for a spherical system [201] . For an equilibrium 
shape that possesses all the rotational degrees of freedom 
of a three-dimensional body (and hence violates spheri¬ 
cal and axial symmetry) such a rotational band on top 
of every intrinsic state involves (2 J + 1) levels with to¬ 
tal angular momentum J. Each of these levels is itself 
(2 J + l)-fold degenerate, corresponding to the different 
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m-substates |20l| . The rotation induced collective en¬ 
hancement may fade out at higher energy, but this will 
happen only far above the energy range around S n , of 
interest here [Hj, |99l, l205| . As usual @, ll94L ll9a j. it 
is assumed in Eq. (l24l) that the excitation energy E x 

is large as compared to the rotational energy h 
and that the spin projections m are normally distributed 
around 0. Obviously the exponential spin cut off factor 
can be neglected in the limit for small J. The reduction 
by a factor 4 is related to the invariance with respect 
to rotations by 180° about each of the three principal 
axeSj valid for the most general quadrupole deformation 
[sl l20l| . If the nuclear body is completely symmetric 
about one axis, i.e. axially deformed, a decrease of the 
level density by \fnj2a = 8 is expected in the limit of 
s mall .7, and a reduction by a factor up to ~ 300 from 
y/nj2a 3 - arises in the comparison between Eq. IMJ) 
and Eq. valid fo r the level density of completely 

spherical nuclei inii mi. The size of these factors in¬ 
dicates that the dependence of the absolute level density 
on the symmetry of the nuclei is appreciable, whereas the 
amount of deformation enters only via the spin cut off a 
and this can be neglected for spins below 1=1. As the de¬ 
formation parameters /3 and 7 are subject to fluctuations 
around their mean values [J a shape sampling equivalent 
to the one described for the axis lengths may be applied 
to the moments of inertia. Such sampling does not af¬ 
fect the symmetry class and its effect on p(E, J 77 ) was 
tested to be negligible by approximate calculations. Low 
energy oscillations in heavy nuclei have been assumed 
to le ad to vibrational bands and to cause an additional 
mi nn collective enhancement of p(E, J 77 ). In our 
scheme extra level density due to quadrupolar modes is 
accounted for by the breaking of axial symmetry. When 
the effect of low energ y qu adrupole vibrations was in¬ 
vestigated as proposed [20 ll ] with strengths derived from 
the isoscalar quadrupole giant resonance (ISGQR) us¬ 
ing hv V ib = \Exsgqr{ 2 + ) and / m b = \Iisgqr{ 2 + )- In 
analogy to the expression for the strength functions as 
compiled in Table HI an enhancement resulted to be far 
below the effect of triaxiality. 


Comparison to data 

Average level distances D(E,J' K ) = l/p(E,J 7T ) are 
available from ensembles of discrete levels with known 
spin populated in nuclear reactions at low E x and at the 
neutron separation energy S n from compound resonances 
in the cross section for neutron capture in the eV and keV 
range. To demonstrate the energy dependence of the level 
density formalism presented here, results for 81 Sr, 113 Cd 
and 235 U are given in Fig. [HU For the plots odd-n nuclei 
with a satisfactory number of levels were selected from 
the data compiled within RIPL-3 ^l|. The experimental 
information on p(Z,N,E — x, J 77 ) stems from counting 
bound levels up to an E x , above which completeness is 
no longer assured, as well as resonances just above S n 


[3TI Il86| . As their spins are known, only two obvious 
assumptions are needed to obtain uj qp (E x ) from the data 
by using Eq. (l24l) : 

1. Parities are equally distributed in the discrete spec¬ 
trum; 

2. a s in the exponential spin cut-off correction is taken 
from systematics [20Clt |. 

This procedure of applying the spin dependent factors to 
the data allows different spins to be shown in the same 
plot. As will be shown in section H the slope below S n 
is of special interest for capture cross sections. It is rea¬ 
sonably well predicted even for 81 Se, although some data 
are missed in absolute height. Extrapolations into do¬ 
mains away from stability need global parameterizations 
valid for an extended region of nuclei. For 132 nuclei 
with A > 70 the average dis tance of s-wave neutron cap¬ 
ture resonances is available |l43j |. As these all have spin 
l/2 + , it is interesting to compare these data to Eq. (1201) . 
which represents the most general, triaxial case. It is 
worthwhile noting that for spin 1/2 the small J limit is 
correct within a few per cent only. Fig. [HU depicts the 
results from using Eqs. fl9l to l23l and l24l in a comparison 
to measurements; obviously many of them lie close to the 
prediction. The shell correction plays an important role 
for the resulting agreement; this is improved by including 
a damping of it as proposed on the basis of an extra 
entropy change with temperature. Agreement over many 
orders of magnitude is reached with the only fit param¬ 
eter a = 0.06 and it is remarkable that even near 208 Pb 
our ’’triaxial” ansatz does quite well. Some alternatives 
proposed earlier may eventually help to obtain an even 
better global fit: The reduction of the nucleon mass to 
an effective in-medium value was suggested [13 , some 
of the expected deviations of nuclear radii from 1.2 • A 1 / 3 
fm are included as calculated if, model calculations of 
moments of inertia show considerable differences to rigid 
rotor values [13, and ambiguities may be seen in the 
choice of temperature (T vs.t). These topics may add 
some to the collective enhancement given by Eq. I!H1) . 
but an important influence on the intrinsic state density 
u> wp (E x ) as well as on p(E x , J) was found to emerge from 
the choice made for 5E, by which the Fermi gas zero is 
fixed with respect to the nucleus’ ground state: Replac¬ 
ing the liquid drop model masses from ref. l208| by the 
one from [203f | increases the level density for actinide nu¬ 
clei by nearly an order of magnitude, wherea s the use 
of a new LDM fit including a curvature term m lias 
the opp osite effect. The recent liquid drop model fit to 
masses [202l | used in this study is based on fitting a vol¬ 
ume and a surface term independently and the resulting 
agreement to resonance spacings are shownjn Fig. [2H1 A 
nearly equally good agre ement is found 1 1241] for a choice 
of LDM parameters [208j | made many years ago. Unfortu¬ 
nately these LDM fits are made neglecting the breaking of 
axial symmetry, as its influence on ground state masses 
was calculated to be very small for most heavy nuclei 
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FIG. 28. (Color online) For the state density to qp (E x ) in the nuclei 81 Sr (a), 113 Cd (b) and 235 U (c) the prediction is shown 
together with respective data from RIPL-3 |3l| obtained using discrete levels (black squares) [HI, 1 143ll as well as neutron 
resonance spacings (below vertical arrow); both are converted into state density u> by inverting Eq. (1241) . A change in slope 
at the phase transition energy E p t is clearly seen (diagonal arrow). The dotted blue lines depict results from Eq. (1191 and l20l) 
using parameters as given in Eqs. (En and EH}. 



FIG. 29. (Color online) Average resonance spacings 
D{S n l/ 2+) in nuclei with 51 < A < 253 as observed in 
neutron capture by spin 0 target nuclei. Data (black data 
symbols) compiled by the RIPL collaboration [3ll . Il43l1 are 
comp ared to the prediction with an effective shell correction 
[2021 ] as presented in section Hxl f dotted curve in magenta). 
The drawn blue line depicts the case when it is damped with 
E x , as proposed previously [ 2041] . 


zMH- (fence more studies on the 4-dependence of SE 
seem indicated; a broad mass range should be regarded, 
as a large effect of the LDM fits to 6E is seen for A ss 60 
and A > 230, but only small differences appear between 
various LDM for A « 80 to 140. The agreement to ex¬ 


perimental level densities from using this LDM fit [202!] 
may well be regarded as an additional indicator for its 
quality, and it especially supports our parameterization. 


X. RADIATIVE NEUTRON CAPTURE 

General remarks 

The radiative capture of fast neutrons by heavier nuclei 
plays an important role in considerations for advanced 
nuclear systems and it is of interest also for the cosmic 
nucleosynthesis. To test the combination of the present 
ansatz on the photon strength to the one for the level 
density - both allowing for a breaking of axial symmetry 
- a comparison on absolute scale of predicted to mea¬ 
sured average radiative widths is demonstrated first. A 
sum over the decay channels to all bound states Jb which 
can be reached from the capture resonances J r by pho¬ 
tons of energy E 1 = E r — Eb, multiplied by their density 
p(Eb, Jb) leads to an effective averaging. The dependence 
of p(Eb, Jb) on Jb ( cf . Eq. (124} ) has the consequence 
that the y-decay for dipole-transitions from J r = 1/2 
to Jb = 3/2 is favored as compared to final levels with 
Jb = 1/2 and a difference between El and Ml arises, 
if there is a parity dependence of the level density. In 
the ansatz discussed here, p(E x ,J) above E pt is deter¬ 
mined by the Fermi gas prescription and independent of 
parity by principle. One may account for differences in 
in the low energy regime caused by a dominant parity 
through a respective estimate of w(0) = u(E x < A 0 ). If 
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respective information is available this may improve the 
predictions. The mean radiative width of capture reso¬ 
nances with J x = l/2 + is the basis for the description 
of radiative capture by 0 + -target nuclei with l n = 0 as 
discussed in the following. 


Average radiative widths 


It was pointed out previously 115] that strength in¬ 
formation can be extracted from capture data directly 
by regarding average photon widths f 7 . These are pro¬ 
portional to the ratio between the level densities at the 
capturing resonances r and at the final states b below S n 
reached by Ej = E r — E b . Eq. (fl4l) shows that the first 
mentioned p(E r ,J r ) are included in fi(E-y). Of course, 
the widths T 7 depend in addition on the photon strength, 
i.e. the low energy IVGDR tail extrapolated fro m ab ove. 
It is known that T 7 does not vary with E r EMM and 
hence it can be approximated for J r = l/2 + by setting 
A = 1 in Eq. (1141) and summing over all final bound levels 
b £ Af „i.e. over A b = [0, S n + E r ]: 


f 7 = E 


be A b 


= f p(E b , J{,)(T 7 (r -a b))dEj 

J A h 


(25) 


r p(E bl j b ) 

IA b P(.Er, Jr) 


h(E^dE y 


Average radiative widths were derived by a resonance 
analysis of neutron data taken just above S n and tabu¬ 
lated [l43j for 115 even-odd nuclei with 51 < A < 253. 
These T 7 allow for a combined test of predictions for 
photon strength and level density slopes, and respective 
data are shown in Fig. [30] Radiative capture into spin-0 
targets through the s-channel only (l n = 0) is consid¬ 
ered, such that the spin depende nce o f the level densities 
p(E r ,J r ), under study recently |200j . can be neglected 
and the spin cut off in Eq. (l24ll is neglected here as the 
collective energy to be added to the intrinsic excitation 
is small as compared to E 1 . Magnetic (Ml) strength 
was treated separately and then the number of levels 
with J n = l/2 + and 3/2 + at low energy enters explic¬ 
itly. They play a minor role for the energy integrated 
photon interaction, as resulted from inserting tabulated 
{HI, [31] values for the matrix elements and recent ex¬ 
perimental data pies eh. It appears that, similar 

to non-nuclear systems, electric dipole modes dominate 
radiative processes. 

In Fig. [30] the agreement on absolute scale between 
prediction and data appears as satisfactory and this can 
be attributed to the combination of TLO to a LDM 
based shell correction. There are some discrepancies 
probably related to the neglect of shell effects (in ad¬ 
dition to SE(Z,A)). In Fig. [301 th e photon strength 
as discussed in sections |V! to lYIIll and the level den¬ 
sity obtained from the parameterization described in sec- 



FIG. 30. (Color online) Average radiative widths as deter¬ 
mined experimentally for nuclei populated in radiative neu¬ 
tron ca pture by spin 0 targets with A > 50 shown as black 
squares [l43ll are compared to calculatio ns. Th e combination 
of the dipole strength (see sections 0 to lVIlJ) with the level 
density I'sectior lTXl) is depicted as drawn line in blue and the 
dotted line in magenta indicates an increase of u>u,p(0) by a 
factor of 3. 


tion llXI are used. Previous work, which only covered lim¬ 
ited numbers of nuclides, used parameters for p(E x ) lo¬ 
cally adjusted to low lying levels and neutron resonances 
[M em em eh Eg. The good overall accordance 
to experimental data reached by us tests the various ap¬ 
proximations applied: A decrease of w,„p(0) by a factor of 
3 modifies f 7 by 20 to 50% when regarding nuclei with 
Am 70 resp.A ss 240; this shows that the agreement 
may be improved by an introduction of local spectro¬ 
scopic information in addition to our parameterization, 
explicitly based on global properties only. The inclusion 
of all the various components of minor strength as listed 
in Table [U increases f 7 by ss 55% for an average over A. 
The largest effect (about 30%) was found for the cou¬ 
pled mode (2 + x 3 - )x-, studied theoretically since 
long. For its strength only data scattered in A and Z 
are available [18] and they were used as guide here. The 
second largest increase of P 7 (by m 10%) and more for 
nuclei with large Q 0 is due to magnetic dipole strength 
(of the scissors mode). A low-energy E2 component of 
strength as discussed in the context of Table [I] was found 
to increase the decay of the l/2+-levels just above S n , as 
populated via s-capture, by less than 7% albeit the open¬ 
ing of E2-decay to 5/2 + levels. The inclusion of pigmy 
strength seems to have a small effect when an average 
over A is regarded: it stays below 10% for each of the 
two components. 
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Averaged cross sections for radiative capture 

Above s=s 1 keV the neutron widths are that large, that 
(r n ) T 7 and the average width ratio in Eq. (15B1) can 
be replaced by F 7 with an accuracy of a few percent. 


This is known from measured neutron strengths as listed 
for RIPL-3 (3ll . ll43| . Approximating transmission coeffi¬ 
cients (from optical potentials) for the compound nucleus 
formation and a full statistical treatment of its decay by 
a schematic ansatz (209| one finds: 


r r 

-L r?.-L ■ 


Tnlb 


(o~(rc,7))r = 2n 2 XlJ2( 2l + \ f )r-p(E r ,J r ) with ( p =r ) r = T 


r 4 - r 

x n w i 7 


As r 7 is contained implicitly in fi(E 7 ) one gets from Eq. (15B1) cr capt (E n ) for E n = E r — S n : 

pEr 

o-ca P t{E n ) = (cr(n, 7 )) r = 27r 2 A 2 ^(2/ + 1) • / f 1 {E 1 )E^p{E b ,J b )dE 1 

i,h 


(26) 


(27) 


In this cross section formula for cr cap t(E n ) the neutron 
energy enters mainly via A„ and the neutron angular 
momentum l. Following a schematic compound nucleus 
model [2091 ]. one arrives at Eq. (1281) for the Maxwellian 
averaged radiative capture cross section (MACS), aver¬ 
aged over many resonances. The case of non-zero target 
spins in radiative capture as well as the inclusion of l > 0, 
direct capture and inelastic scattering are not treated 
here; respective numerical tests indicate small changes 
only to our conclusions. A prediction of MACS b ased o n 
the photon strength presented in sections [v3 to IyH is 
then possible and calculated by: 


(a-(n,7)) fcT 


2 J™ a capt (E n )E n e~ E ^ kT dE, 
f 0 °° E n e~ En ! kT dE n 


Here cross sections given by Eqs. [TB] and [57] are folded 
with a Maxwellian distribution of neutron energies. As 
pointed out [ 21Cli |. the same can be done with experi¬ 
mental data. In view of the fact that D T r > T r7 
the Maxwellian averages around 30 keV are formed inco¬ 
herently with neglect of Porter-Thomas fluctuations. By 
only regarding the radiative capture by spin-zero targets 
effects related to ambiguities of spin cut off parameters 
and angular momentum coupling are suppressed. It is ob¬ 
vious from Fig. 1311 that the TLO-parameterization with 
the proposed ansatz for p(A, J ™, E x ) work well: For the 
discussed range of A the overall agreement on absolute 
scale covering more than three decades is remarkable. 
The discrepancy observed in the region of A > 230 and 
A < 70 may be related to a false prediction of minor 
strength, or of shell correction influencing level densities, 
and details of the shell structure in these nuclei may be 
important for both. Interesting aspects become obvious 
from the regard of Figs. I5B1 and [311 together: Near closed 
shells like at 208 Pb P 7 is large whereas a minimum ap¬ 
pears for (<j(n, 'y))kT- The average radiative widths di¬ 
rectly rise with fi(E^) and thus E 0 — S n and vary only 
with the slope of p(E Xl J) in the range from E b to E r . 
In contrast, capture cross sections directly depend on the 


level density in the final nucleus p(E b ,J b ), and this re¬ 
sults in a deep minimum at shell closure ( cf \26\\27\ i . From 
this simplified study the main topic of our approach - the 
importance of triaxiality in heavy nuclei - is obvious and 
not influenced by the rather small discrepancies seen in 
Figs. [30] and EU 



FIG. 31. (Colo r on line) Maxwellian averages of experimental 
cross sections [2 111 . l212tl for radiative neutron capture into 
even nuclei with J = 0 and 50 < A < 250 shown as black 
data symbols for kTAGB = 30 keV. They are plotted vs. Acn 
in comparison to calculations based on Eq. (1281) with TLO 
(dotted line in magenta) and including the minor components 
listed in Table Q] (full blue line). The level density prediction 
is described in section HXl a shell effect damping is included. 

It was worked out earlier [83| as well as in section Ivnl 
that the impact of photon strength on radiative neu¬ 
tron capture cross sections is peaking in the region below 
Ery = 6 MeV. Below, the factor i? 2A+1 reduces the tran¬ 
sition rates and above the density of levels to be reached 
becomes small. Thus the scissors mode has to be re¬ 
garded for Ml-transitions in deformed nuclei, and the 
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low energy El-strength in the (2 + x 3“)i_ and the lower 
pigmy modes have a strong effect in all nuclei. As de¬ 
picted in Fig. [31] by the dashed curve, the cross section 
is dominated by the IVGDR extension to lo w ene rgy, es¬ 
pecially for A < 150. Recent work [32j, [zl, Il36( | on the 
influence of minor strength components has indicated a 
more significant impact in comparison to the IVGDR tail, 
which was quantified in a different way; it was shown in 
section IVII1 that single or double Lorentzians result in 
erroneous estimates of the corresponding El-strength in 
the tail. This is especially so if a dependence of P 7 on 
E 7 is imposed [ 8 J, |95[ to seemingly improve the fit of the 
IVGDR. Namely the ganrma-energy dependence of the 
resonance widths imposed often plays an important role. 
In the TLO approach the sole variation of Ti with the 
pole energies its global adjustment and the strict im¬ 
plementation of the TRK sum rule avoid irregularities in 
the A-dependence of the dipole strength. We have shown 
that this is appropriate in the valley of stability, but for 
very neutron rich nuclei with their small S n , experimental 
tests may became possible eventually in newly available 
radioactive beam facilities. 


Fast neutron induced transmutation 

Neutron capture processes play an important role in 
the study of cosmic nucleosynthesis, especially for sce¬ 
narios with such high neutron flux, that neutron cap¬ 
ture processes can be the cause for a production of nu¬ 
clides beyond Fe [H, |24[ . The radiative capture of fast 
neutrons by heavier nuclei is of interest also for consid¬ 
erations for advanced nuclear systems [|(1 | 2 l| and de¬ 
vices aiming for the transmutation of radioactive nuclear 
waste. One for many even the strongest - of the ar¬ 
guments against a long-lasting commitment to nuclear 
power as energy supply is the difficult permanent dis¬ 
posal of the long lived radioactive waste produced in nu¬ 
clear reactors. Significant efforts are thus being made 
worldwide directed towards the minimization, manage¬ 
ment, and disposal of highly radioactive nuclear waste. 
Thus the partitioning of nuclear waste and transmuta¬ 
tion of long-lived isotopes to nuclides with shorter life¬ 
time are being investigated worldwide. Different schemes 
have been proposed that may reduce the radioactivity 
and radio-toxicity of the spent fuel after burn up. The 
studies towards the choice of the best options make ex¬ 
tensive use of simulation methods in order to predict the 
system behaviour in a great variety of possible configura¬ 
tions and running conditions [22] . A fundamental prereq¬ 
uisite for these computations is the availability of reliable 
cross section data. These are needed for processes and 
operating parameters that are significantly different from 
those of current systems mostly operating with thermal 
neutrons. A good understanding of the radiative capture 
of fast neutrons is an important part for any theoreti¬ 
cal attempt in this direction, and the agreement shown 
in Fig. [31] appears promising. But a much more severe 


problem with nuclear reactor waste is posed by the pres¬ 
ence of nuclides beyond Uranium produced in the high 
flux of slow neutrons captured by the 238 U usually form¬ 
ing a large fraction of the fuel. By multiple neutron cap¬ 
ture and subsequent /3-decay an appreciable amount of 
long lived nuclides are formed. As several of these “mi¬ 
nor actinides” are a-emitters the risk of their incorpora¬ 
tion increases their high radiotoxicity. Attempts to ex¬ 
pand the TLO-based calculations of (n,y) cross sections 
to actinide nuclei are hampered by serious discrepancies 
between IVGDR data observed in different laboratories, 
some of which are obvious from figures shown here. Neu¬ 
tron capture by actinide nuclei is of great importance for 
the transmutation of nuclear waste and a dedicated in¬ 
vestigation of neutron capture cross sections for Th, U 
and heavier nuclei, for which data were compiled [ 212 | , is 
useful. The fact, that TLO results in a good description 
of the data for 239 Pu (c/. Fig. [25j and the good match 
observed in Fig. [31] leads to consider the TLO approach 
as encouraging to serve as guideline for predictions con¬ 
cerning the radiative capture of fast neutrons also for 
actinide nuclei. The approximations made to arrive at 
Eqs. [25] and [27] work well in the range of E n ss 30keV 
(see Fig. EU, but the coupling to other channels like in¬ 
elastic scattering has to be included at higher energies. 
The effect of further contributions was shown recently in 
detail for 238 U 621 and for other actinides |213j . where 
also the importance of the scissors mode was pointed out. 
In view of the average overall agreement for more than 
150 nuclides an inclusion of the TLO dipole strength and 
the “triaxial” le_vel density as an option in TALYS jl86l | 
or EMPIRE [214l | appear highly desirable. These two 
codes now use the “old” Lorentzian descriptions [31], E3] 
of the IVGDR shapes with only one or two poles and 
locally adjusted IVGDR parameters. 

The extension of our ansatz to targets with nonzero 
spin needs additional caution as the transition to the ob¬ 
server’s system in Eq. (ED assumes the rotational energy 
to increase with spin. Data for resonance spacings near 
S n as observed with neutron capture by odd nuclei in¬ 
volve an elevated spin of usually single particle nature. 
It is thus disputable to simply apply Eq. (1241) to higher 
spin, even if one of the several propositions [ 200 | for the 
spin cut off parameter a is inserted. A related ambiguity 
emerges, when one of these is inserted into a factoriza¬ 
tion for the spin dependence: p(E x ,J) = S a (J)p(E x ). 
The ansatz presented here for J= 1/2 avoids this prob¬ 
lem and further study is needed to extend it to higher 
spin. 


XI. CONCLUSIONS 

Various spectroscopic information presented over the 
years ai, min I47l ] indicated triaxiality for a num¬ 
ber of heavy nuclei. Admission of the breaking of axial 
symmetry, albeit often weak, clearly improves a global 
description of Giant Dipole Resonance (IVGDR) shapes 
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by a triple Lorentzian (TLO), introduced recently fl3l. |72| 
and discussed in detail by this paper. The three parts 
add up to the TRK sum rule |53l [54J, when theoretical 
predictions for the A-dependence of pole energies from 
droplet model |70|_and spreading widths based on one- 
body dissipation [74| are used. Two additional effects 
- hitherto not emphasized as such - indicate a break¬ 
ing of axial symmetry in nearly all of the heavy nuclei: 
With one global parameter the scheme proposed by us 
reproduces observations for level densities in nuclei with 
A > 50 and J = 1/2, when (a) the Fermi gas prescription 
is only used above a phase transition at the critical tem¬ 
perature, (b) a pairing condensation energy is included 
in the backshift and (c) the collective enhancement due 
to symmetry breaking by triaxiality is included. Only 
a modest modification of the level density parameter a 
from its nuclear matter value is needed to fit resonance 
spacing data; this is a surprisingly small surface term , 
the only free parameter of our level density prescription. 
Again only one global free parameter suffices to fit to the 
shape of the IVGDR peak by a triple Lorentzian pho¬ 
ton strength (TLO), somewhat improved as compared to 
the original proposal El lil , well in accord to the TRK 
sum rule. It also predicts its low energy tail - without 
other modification than the addition of minor strength 
components - to match respective strength data. This 
consideration of broken spherical and axial symmetry - 
for low excitation in accord to CHFB calculations E] and 
even more so for increasing excitation energy - allows for 
a combination of both parameterizations. Exact defor¬ 
mation parameters are unimportant for the tail of the 
El-resonance as well as for the density of low spin states 
occurring in neutron capture by even targets, where nei¬ 
ther spin cut off nor moments of inertia are involved. 
Thus a combination of the ansatz for photon strength to 
the novel approximation for level densities leads to a pre¬ 
diction on absolute scale for neutron capture in the range 
of unresolved resonances - including Maxwellian average 
cross sections compiled recently [211] for ( E n ) = 30keV. 

Four points are stressed: 

1. Ad hoc assumptions on shapes and collective en¬ 
hancement in previous calculations (3l) of com¬ 
pound nuclear rates are replaced by a direct ac¬ 
count for broken symmetries. 

2. Doing so, the prediction for level densities is im¬ 
proved such, that only one small parameter suffices 
to obtain agreement to data on absolute scale. 

3. The low energy tail derived from the global triple 
Lorentzian (TLO) fit to IVGDRs is not modified 
by an extra energy dependence. 

4. The predictions for average radiative widths and 
Maxwellian cross sections are sensitive to one in¬ 
gredient to the calculations not controlled globally: 
The state density co qp ( 0) near the ground state of 
the final nucleus. 


5. A similar statement has to be made concerning the 
low energy electric dipole strength induced by cou¬ 
pling nuclear quadrupole and octupole modes. For 
these a glo bal prediction appears impossible and for 
Figures |7j to [3l| values were used, which were de¬ 
rived as upper limits to experimentally determined 
data. 

6 . The possible increase of n-capture yields by some 
amount was discussed [3 to be cause d by addi¬ 
tional photon strength, predicted EH to appear 
below ft: OAEivgr ■ Such pigmy strength may 
eventually be related to a vibration of excess neu¬ 
trons against a core. Another excess above the 
Lorentzian tail near S n seem s to be related to a 
vortical proton motion [l79j |. Our study finds a 
rather weak influence of such strength components 
in comparison to the one mentioned under 5. 

In conclusion: For more than 100 spin-0 target nuclei 
with A > 50 resonance spacing data and average cap¬ 
ture cross sections are well described by a global fit with 
only a small number of adjusted quantities, which turn 
out to be A-independent. This successful ansatz can be 
used as basis for detailed studies in specific regions of the 
nuclide chart and this promises reliable predictions for 
other compound nuclear reactions besides radiative cap¬ 
ture. The breaking of axial symmetry in excited heavy 
nuclei is demonstrated here on the basis of experimen¬ 
tal data, most of which were published previously. Fur¬ 
ther studies on broken axial symmetry in addition to ex¬ 
isting theoretical hints to the im p orta nce of triaxiality 

El SB S3 m m 0 0, m El El mu seem neces¬ 
sary, as nuclear theory studies often still prefer to assume 
axial symmetry ad hoc. We actually question that an ex¬ 
perimental proof for such symmetry exists, and probably 
all heavy nuclei do not strictly obey it, in apparent resem¬ 
blance to what was found for crystalline matter by Jahn 
and Teller in 1937 [216j . Such an analogy is suggested 
by the data analysis presented here, showing that the 
global ansatz derived on the basis of broken axial sym¬ 
metry combined to the dipole sum rule allows remarkably 
good predictions for radiative neutron capture in general 
- with very few free fit parameters. 
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